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Final  Report  for  Developing  Equalization  Techniques  for  use  with 
Bandwidth  Efficient  Modulation  Methods 

Subtitle:  Data  Compression  Techniques  to  Reduce  Bandwidth 


Abstract 


Bandwidth  is  a  precious  commodity.  In  order  to  make  the  best  use  of  what  is  available, 
better  modulation  schemes  need  to  be  developed,  or  less  data  needs  to  be  sent.  This 
paper  will  investigate  the  option  of  sending  less  data  via  data  compression.  The  structure 
and  the  entropy  of  the  data  determine  how  much  lossless  compression  can  be  obtained  for 
a  given  set  of  data.  This  paper  shows  the  data  structure  and  entropy  for  several  actual 
telemetry  data  sets  and  the  resulting  lossless  compression  obtainable  using  data 
compression  techniques. 


Final  Report  for  Developing  Equalization  Techniques  for  use  with 
Bandwidth  Efficient  Modulation  Methods 

Subtitle:  Data  Compression  Techniques  to  Reduce  Bandwidth 

Project  Summary; 

Objectives: 

1 .  Analyze  the  performance  of  the  32  QAM  in  multipath  channels. 

2.  Develop  adaptive  filtering  equalization  for  32  QAM  in  multipath  channels. 

3.  Analyze  the  performance  of  pulse  shaped  8  PSK  in  multipath  channels. 

4.  Develop  adaptive  filtering  equalization  for  8  PSK  in  multipath  channels. 

5.  Investigate  the  use  of  interleaving  for  the  32  QAM  and  8  PSK  modulation  techniques 
in  a  multipath  environment. 

6.  Develop  reports  detailing  work  and  recommendations. 

7.  Perform  statistical  analysis  on  the  telemetry  data  to  determine  possible  amounts  of 
compression. 

8.  Apply  arithmetic,  Ziv-Lempel,  Huffman,  and  Rice  data  compression  techniques  to  the 
data. 

9.  Test  data  compression  techniques  in  a  channel  simulation  to  test  the  robustness  of  the 
code. 

10.  Adapt  the  Rice  algorithm  to  better  handle  telemetry  data. 

11.  Complete  Final  report. 

Accomplishments: 

Objectives  1  through  6  have  not  been  met  and  were  changed.  To  test  the  32  QAM  and  8 
PSK  in  the  multipath  channels  required  the  appropriate  channel  model  be  available. 
Brigham  Young  University  (BYU)  was  developing  this  model  when  this  project  started, 
and^the  model  was  not  available.  During  the  progress  of  this  project,  the  objectives 
changed.  Discussion  with  ARTM  led  to  a  different  area  of  concentration.  Work  by 
others  was  proving  successful  with  FQPSK,  so  this  work  on  modulation  was  changed  to 
an  emphasis  on  data  compression.  Reducing  the  eimount  of  data  to  be  sent  can  reduce  the 
bandwidth  and  thus  accomplish  the  same  results  as  the  8  PSK,  and  32  QAM. 

Objectives  7  through  9  were  completed.  Data  compression  results  indicate  that  data 
compression  is  viable,  channel  simulations  do  indeed  indicate  that  data  compression 
alone  would  be  unfeasible,  but  adding  channel  coding  to  protect  the  data  again  makes  this 
data  compression  with  channel  coding  a  viable  option. 

Objective  10  is  still  in  process.  The  work  was  not  completed  due  to  the  extensive  amount 
of  time  it  would  take  to  complete  the  task,  and  student  employment  problems.  Current 
results  do  indicate  that  an  adaptation  of  the  Rice  algorithm  would  work  well  for  telemetry 
data. 
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Personnel  Support: 


1 .  One  Faculty  member  (the  PI) 

2.  One  technician  (for  help  in  setting  up  the  computers  in  the  lab) 

3.  Five  undergraduate  students  (total  at  various  times,  no  student  worked  more 
than  part-time.  Two  students  worked  for  lyear,  and  2  other  students  worked 
for  6  months.  The  last  student  was  supported  by  an  on  campus  research  grant 
for  2.5  months)  The  students  were: 

Jonathan  McMillan  -  worked  on  learning  simulations  and  getting  simulink  to 
work;  attended  1998  ITC  conference. 

James  Woods  -  wrote  programs  for  statistical  analysis,  attended  1998  ITC 
conference. 

Shakti  Davis  -  worked  on  the  error  analysis,  and  channel  simulations, 
attended  1999  ITC  conference 

Kendall  Mauldin  -  began  looking  into  the  Rice  algorithm  and  made 
recommendations  to  modify  it,  attended  1999  ITC  conference. 

Enrique  Sanchez  -  started  modifications  of  Rice  algorithm  but  was  unable  to 
complete  the  work. 

Publications; 


Paper  at  the  International  Telemetering  Conference: 

“Data  Compression  Statistics  and  Implications”,  ITC  October  1999 

Interactions/T  ransitions; 

1 .  Participated  in  Advanced  Range  Telemetry  (ARTM)  group  meeting  in  the 
spring  (1998). 

2.  Consultative/advisory  functions:  Communicate  with  BYU  regarding  our 
related  projects 

3.  Participated  in  Advanced  Range  Telemetry  (ARTM)  group  meeting  at  the  ITC 
(1998,  1999) 

4.  Participated  in  Advanced  Range  Telemetry  (ARTM)  communications  via 
email  and  phone.  New  error  sets  were  to  be  sent  in  January  (1999)  -  none 
received  as  of  this  date. 

New  Discoveries,  inventions,  or  patent  disclosures:  None 

Honors/Awards:  None 
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Introduction 


Bandwidth  is  a  precious  commodity.  It  has  become  imperative  that  bandwidth  is  used 
efficiently.  To  this  end,  data  needs  to  be  reduced,  or  modulation  techniques  need  to  be 
used  t  minimize  the  bandwidth.  Data  compression  (coding)  is  one  technique  to  reduce 
the  bandwidth  needed  to  transmit  data.  The  bit  rate  =  R  =  (#  of  data  bits  sent)/  sec.  The 
bandwidths  for  BPSK  are  (6): 

3-db  bandwidth  0.88  R 

null  to  null  bandwidth  1 .0  R 

As  can  be  seen  if  R  decreases,  the  so  to  does  the  bandwidth.  The  goals  of  data  coding  are 
to: 

•  Reduce  the  amount  of  data 

•  Find  codes  that  take  into  account  redundancy,  structure,  and  patterns 

•  Break  up  data  into  random  groups 

Entropy  is  a  measure  of  the  information  content  of  the  data,  a  measure  of  the  uncertainty 
or  randomness  of  the  data,  and  can  be  used  to  indicate  how  rhany  bits  would  be  needed  to 
encode  the  data  if  the  data  are  independent. 

If  X  is  a  random  variable,  process,  or  data  to  be  coded  (usually  X  is  discrete,  but  doesn't 
have  to  be);  which  takes  on  the  values  Xi,  X2,  ...,  Xn;  and  pi  is  the  probability  that  Xi 
occurs,  then  the  entropy 

H(X)  =  -Y.p,\og(p,) 

all  i 

If  the  log  is  base  two,  then  the  units  will  be  bits  and  will  indicate  the  number  of  bits 
needed  to  code  the  message  without  loss  of  information.  In  order  to  calculate  this 
quantity,  the  all  of  the  marginal  and  joint  probabilities  have  to  be  known.  This  is  not 
usually  possible.  Most  of  the  time,  the  entropy  is  estimated  by  the  first  order 
probabilities  (that  is  the  probabilities  of  the  individual  symbols).  A  second  order 
probability  is  based  on  grouping  two  symbols  at  a  time;  third  order  is  three  symbols  at  a 
time.  etc.  The  true  entropy  is  based  on  the  nth  order  probabilities  as  n  approaches  the 
size  of  the  data  file  to  be  transmitted. 

Independence  of  data  means  that  the  probability  of  (X1X2 )  is  the  product  of  the 
individual  probabilities. 

Example: 

Suppose  that  the  original  data  is  given  by: 

0  1  22  23  444  5  66  7  8  8  8 

To  code  this  directly  in  binary  would  take  4  bits  per  symbol.  The  amount  of  data 
compression  can  be  found  by  using  entropy  and  assuming  independence  of  the  data. 

To  find  the  first  order  entropy,  the  probability  of  each  symbol  would  need  to  be  found. 
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Symbol 

Freq. 

Pi  =  freq/(total  # 
symbols) 

-Pi  log  Pi 

0 

1 

0.0625 

0.25 

1 

1 

0.25 

2 

3 

0.453 

n 

:> 

1 

0.0625 

0.25 

4 

3 

0.1875 

0.453 

5 

1 

0.0625 

0.25 

6 

2 

0.375 

7 

1 

0.25 

8 

3 

0.1875 

0.453 

overall 

16 

1.0 

Data  entropy=2.983 

This  gives  a  reduction  of  about  25%.  The  %  compression  can  be  found  by 

compressed  file  size  (average  coded  symbol  size) 

compression  ratio  =  1 - : - — — : — ; - — - ; — — — ; - 

orignal  file  size  (uncoded  symbol  size) 

If  the  data  are  independent  and  if  the  above  is  an  accurate  representation  of  the  entropy,, 
then  the  best  compression  that  would  be  possible  is  2.983  bits/symbol. 

The  following  table  shows  the  results  obtained  if  the  data  are  grouped  by  two  symbols  at 
a  time: 


Two  bits 

New  symbol 

Freq 

Prob 

-Pi  log  Pj 

01 

0 

1 

0.125 

0.375 

22 

1 

1 

0.125 

0.375 

23 

2 

1 

0.125 

0.375 

44 

J 

1 

0.375 

45 

4 

1 

0.125 

0.375 

66 

5 

1 

0.125 

0.375 

78 

6 

1 

0.125 

0.375 

88 

7 

1 

0.125 

0.375 

totals 

8 

1.0 

3.0 

This  gives  an  entropy  of  3.0  (still  a  decrease  from  the  original,  but  an  apparent  increase 
from  the  first  order  entropy).  This  really  isn’t  an  increase.  This  entropy  is  for  two 
symbols  at  a  time  and  the  compression  ratio  will  be  l-(3x8)/(4xl6)  =  63%. 

These  entropies  were  calculated  assuming  independence  of  the  data.  Suppose  there  is  a 
relationship  between  the  data.  One  way  to  examine  this  is  to  look  at  the  differences 
between  neighboring  data.  To  start,  leave  the  first  data  symbol  alone,  and  then  subtract 
the  second  from  the  first,  then  the  third  from  the  second,  etc.  This  gives; 

Original  data;  0122234445667888 
Differenced;  0110011001101100 
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In  this  case,  there  are  two  possible  symbols  ( 0  and  1),  each  occurring  8  times.  This  will 
give  a  probability  of  Vi  for  each  symbol,  and  an  entropy  of  1 .  We  could  now  code  the 
original  data  with  16  bits  of  data  where  the  original  would  take  4  times  that  many  bits. 
This  would  give  a  compression  ratio  of  75%  over  the  original  data.  Since  we  get 
compression  by  differencing  the  data,  this  indicates  that  there  is  correlation  between  the 
data  and  that  they  are  not  independent.  Suppose  now  that  the  differenced  data  are 
grouped  two  at  a  time  as  follows: 

01  10  01  10  01  10  11  00 


Analysis  of  this  data  gives 


Two  bits 

New  symbol 

Freq 

Prob 

-Pj  log  Pj 

01 

0 

3 

0.375 

0.531 

10 

1 

3 

0.375 

0.531 

11 

2 

1 

0.125 

0.375 

00 

3  n 

1  1 

0.125 

0.375 

totals 

8 

1.0 

1.811 

The  compression  ratio  for  this  grouping  is:  1-  (8x1.81  l)/(4xl6)  =  77% 

Notice  how  the  compression  ratio  continues  to  increase,  meaning  that  the  file  size  gets 
smaller.  Continuing  to  group  one  more  time  -  grouping  by  fours  gives: 

0110  0110  0110  1100 


This  has  just  two  symbols.  0110  occurs  three  times,  and  1 1 00  occurs  only  once.  This 
will  give: 


Two  bits 

New  symbol 

Freq 

Prob 

-Pi  log  P, 

0110 

0 

J 

0.75 

0.311 

1100 

1 

1 

0.25 

0.5 

totals 

4 

1.0 

0.811 

In  this  case,  with  just  two  symbols  to  code,  the  compression  ratio  becomes: 

1-  (4x0.811)/(4xl6)  =  95% 

The  above  is  just  an  example  to  demonstrate  how  compression  can  occur.  The 
compression  ratio  is  just  a  prediction  of  the  actual  compression.  In  order  to  achieve  this 
predicted  compression,  an  appropriate  compression  algorithm  would  need  to  be 
developed.  But  this  can  act  as  a  guide  for  how  well  the  compression  could  do.  If  all  of 
the  probabilities  were  known  for  the  data,  then  entropy  would  give  a  lower  bound  on  the 
amount  of  compression  to  be  expected.  As  higher  order  probabilities  are  used,  the 
predicted  compression  will  become  closer  to  the  entropy  bound. 

The  above  ideas  of  data  compression  and  entropy  will  now  be  applied  to  Telemetry  data. 
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Test  Data 


Thirteen  actual  telemetry  data  files  were  obtained  from  Vem  Diekman  at  TYBRIN  (1) 
through  the  Advanced  Range  Telemetry  (ARTM)  group.  The  data  were  then  analyzed  to 
find  the  4  bit,  8  bit,  and  12  bit  entropies.  The  entropy  for  a  set  of  data  is  given  by: 

\  m  m  m 

H{S)  =  -  lim-  S  I ...  E {P{X,  =  ;•  ,^2  =  =  /„)  • 

=1/2=1  /„=i 

=  /j,  J^2  ~  ~  ^"«)} 

(1) 


Where  H  (S)  is  the  entropy  of  the  source  (the  data) 


Xj  are  the  elements  of  all  possible  sequences  of  the  data.  The  X,  can  be  bits, 
symbols,  or  words. 

n  =  the  length  of  the  sequence 
P  is  the  probability 

m  =  the  size  of  the  alphabet  used  for  coding;  for  binary  data,  m=2,  for  words  of  4 
bits  in  length, 

m=16,  etc. 

To  find  the  4-bit  entropy,  all  possible  combinations  of  four  elements  from  the  binary 
alphabet  were  found.  The  frequency  of  occurrence  of  each  of  these  4  bit  words  was 
found.  This  frequency  was  used  as  an  estimate  of  the  probability  of  occurrence  for  each 
of  these  words.  From  these  probabilities,  the  estimate  for  the  4  bit  entropy  was  found 
using  equation  1  where  n=4  and  m=2.  For  the  8-bit  entropy.  256  different  words  are 
possible,  and  for  the  12-bit  entropy,  4096  words  are  poss’bie.  It  can  be  seen  that  letting  n 
approach  infinity  will  quickly  become  impractical.  The  formats  for  each  data  set  varied. 
Some  of  the  data  were  coded  into  words  of  10  bits,  some  12  bits,  etc.  Plots  and  statistics 
for  each  of  the  data  sets  are  in  the  appendix.  A  summary  of  the  data  statistics  is  given 
below: 


Data  Set  Entropies  and  Word 

Size 

Data  Set 

Data  word  size 
In  bits 

4  bit 
Entropy 

8  bit 
Entropy 

12  bit 
Entropy 

#bits/file 

SDSOOl 

12 

3.9 

7.5 

10.3 

139536000 

SDS002 

12 

3.3 

6.0 

7.5 

94239360 

SDS003 

12 

3.8 

7.0 

9.0 

157352160 

SDS004 

10 

3.4 

5.9 

7.3  1 

496743000 

SDS005 

10 

3.6 

6.6  n 

8.6 

1234710000 

SDS006 

10 

3.7 

5.9 

7.6 

890556000 

SDS007 

10 

2.9 

5.9 

7.5 

2347503000 

SDS008 

12 

2.6 

4.2 

4.7 

1105437696 
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Data  Set 

Entropies  and  Wore 

Size 

Data  Set 

Data  word  size 
In  bits 

4  bit 
Entropy 

8  bit 
Entropy 

12  bit 
Entropy 

#bits/file 

SDS009 

12 

3.2 

5.8 

6.9 

1102886400 

SDSOlO 

12 

2.7 

4.3 

4.8 

1100703744 

SDSOll 

12 

2.8 

4.9 

5.9 

568813056 

SDS012 

16 

2.5 

3.7 

5.2 

2193232356 

SDS013 

24 

2.8 

3.7 

4.6 

3443281920 

To  determine  the  amount  of  compression  that  should  be  possible,  the  entropy  per  number 
of  bits  must  be  used.  Hence  %  compression  for  k  bit  entropies  =  (l-(k  bit  entropy  / 
k))*  1 00.  The  amount  of  possible  compression  for  each  data  set  can  then  be  found.  If  the 
k  bit  entropy  is  equal  to  the  entropy,  then  the  predicted  compression  will  be  a  lower 
bound  for  all  data  compression  techniques.  A  plot  of  the  predicted  compression  from,  the 
calculated  entropies  is  given  in  Figure  1 .  It  can  be  seen  that  the  compression  increases  as 
the  number  of  bits  per  symbol  is  increased.  This  indicates  that  the  entropies  do  not  yet 
equal  the  entropy  for  the  data  sets.  Consequently,  compression  techniques  should  be  able 
to  achieve  values  better  than  these  predicted  values. 

The  Huffman,  an  Adaptive  Huffman,  Arithmetic,  and  Lempel-Ziv  compression 
algorithms  were  applied  using  programs  from  Mark  Nelson’s  text  [2].  The  Huffman 
codes  and  Arithmetic  codes  use  an  8-bit  word  length  in  the  code.  Two  variations  of  the  . 
Lempel-Ziv  algorithm  were  used.  The  LZSS  uses  a  pair  of  values  to  indicate  the  location 
of  the  match  and  the  length  of  the  match  in  the  dictionary.  The  LZW  algorithm  involves 
only  sending  one  element  instead  of  a  pair  of  elements,  and  using  a  start  up  alphabet  in 
the  dictionary  consisting  of  all  the  letters  of  the  source  alphabet.  The  percent 
compression  for  each  technique  along  with  the  predicted  compression  is  given  in  the 
following  table. 


Data  Compression  Resu 

.ts 

Data 

Set 

Huffman 

Adapt. 

Huffman 

LZSS 

LZW 

Arith¬ 

metic 

%  4  bit 
Predicted 
compress 

%  8  bit 
Predicted 
compress 

%  12  bit 
predicted 
compress 

SDSOOl 

6.3 

6.9 

15.5 

-9.9 

6.6 

2.5 

6.3 

14.2 

SDS002 

24.1 

25.5 

61.1 

33.8 

24.7 

17.5 

25.0 

37.5 

SDS003 

11.5 

15.5 

39.7 

-7.3 

MnsM 

5.0 

12.5 

25.0 

SDS004 

25.1 

27.7 

50.2 

11.5 

■aais^ 

15.0 

26.3 

39.2 

SDS005 

18.5 

19.0 

MM 

18.5 

18.7 

10.0 

17.5 

28.3 

SDS006 

25.1 

25.9 

BfiM 

32.5 

25.4 

7.5 

26.3 

36.7 

SDS007 

25.9 

26.5 

41.9 

28.2 

mm 

m^BM 

26.3 

37.5 

46 

48.3 

BiSM 

47.5 

60.8 

MiMiIifel 

27.5 

29.4 

64.7 

26.3 

27.7 

27.5 

42.5 

SDSOlO 

45.2 

47.5 

68.9 

52.7 

45.9 

32.5 

46.3 

60.0 

SDSOll 

37.6 

39.7 

46.0 

38.1 

30.0 

38.8 

50.8 

SDS012 

51.9 

51.9 

51.9 

37.5 

53.8 

56.7 

SDS013 

48.6 

53.7 

65.5 

64.8 

51.4 

30.0 

53.8 

61.7 
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A  negative  with  the  compression  value  indicates  that  the  file  was  expanded  instead  of 
compressed.  It  can  be  observed  that  each  data  file  compressed  differently.  Figure  2 
contains  the  plot  of  the  results  of  the  compression  techniques.  It  can  be  seen  from  Figure 
2  that  certain  files  will  compress  very  well.  A  60%  compression  would  mean  that  the  file 
would  take  up  less  than  half  its  original  size.  In  all  cases,  there  is  at  least  one  technique 
that  provides  compression  of  10%  or  more.  With  the  demands  on  spectra,  even  this  little 
gain  can  be  worth  something.  Since  the  Huffman  and  Arithmetic  codes  that  were  tried 
work  with  8  bit  word  sizes,  if  the  compression  obtained  by  these  techniques  is  compared 
with  the  8-bit  entropy,  we  see  a  very  close  match.  All  but  2  of  the  results  are  within  1 0% 
or  less  of  the  8  bit  entropy.  Since  the  LZ  algorithms  achieve  larger  compression  than  this 
indicates,  the  8-bit  word  size  is  not  the  best  choice  and  that  this  8-bit  entropy  is  not  the 
true  entropy.  Also,  since  the  k  bit  entropies  continue  to  increase  with  k,  this  also 
indicates  that  the  actual  entropy  has  yet  to  be  found.  Hence  even  larger  compressions  can 
be  expected. 

The  Rice  data  compression  algorithm  was  also  tried.  The  results  were  discouraging. 
Instead  of  any  compression,  files  were  either  expanded  or  were  unable  to  decode  to 
regain  the  original  file.  The  main  problem  with  the  Rice  algorithm  is  that  is  was  designed 
specifically  for  images.  The  next  thing  tried  was  to  modify  the  Rice  code  witten  by 
UNM  to  run  with  variable  word  length,  and  variable  number  of  words  per  frame.  This 
was  found  to  be  much  more  work  because  of  the  following  problems: 

1 .  The  code  is  written  and  optimized  to  deal  with  image  files,  not  telemetry  data.  This 
causes  several  “domino-effect”  problems. 

2.  A  “long”  (32-bit  variable)  is  used  to  keep  track  of  the  number  of  pixels  (bits  in  our 
case)  in  a  file.  The  largest  value  possible  for  a  “long”  is  2,147,483,647.  So  if  a  file 
size  ever  exceeds  2,147,483,647  bits  (268,435,456  bytes  =  approx.  256  Megabytes), 
we  would  hit  this  limitation  -  one  telemetry  file  exceeds  this  limit. 

3.  ARTM  telemetry  is  organized  in  frames  with  length  of  200  words  each.  The  UNM 
code  (without  modification)  doesn’t  support  this  frame  size  or  structure. 

4.  As  said  in  a  comment  on  line  3734  of  rice.c  (UNM  code)  the  hardware  (and  thus  the 
C  code  that  simulates  the  chip)  is  limited  to  64  blocks  per  scanline. 

5.  There  also  seems  to  be  a  magic  number  of  64  when  it  comes  to  the  “pixels  per  block” 
command-line  option  of  the  compression  program.  When  this  option  (-j  <integer 
number  here>  )  is  set  to  64  or  less,  there  is  no  problem  compressing  or 
UNCOMPRESSING  the  file.  When  it  is  set  to  >64  there  are  no  problems 
compressing  the  file,  by  when  the  file  is  UNCOMPRESSED,  the  uncompression  part 
of  the  program  gives  an  error  message  and  halts.  Also,  the  UNM  code  truncates  the 
data  frame  if  it  is  at  the  end  of  the  file  and  its  length  is  less  than  a  complete  scan  line. 

Noting  these  problems,  the  task  of  re-writing  the  code  was  undertaken,  but  time  ran  out 
before  the  task  was  completed.  In  an  effort  to  decide  whether  this  work  would  be  worth 
the  effort,  other  work  was  investigated. 
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Figure  1.  Predicted  Compression  from 
Data  Statistics 


-^SDSOOI 
•  -  SDS002 
SDS003 
••  -  SDS004 
-^SDSOOS 
SDS006 
-^SDSOOT 
SDS008 
-4— SDSOOSi 
SDS010 
SDS011 
SDS012 
SDS013 


IT. 


•  SOS002 
•-SOS003 


♦•-SDS005 

^•SDSooe 


••SDS008: 

-sosocsj 

-SDS010 

••SDS011 

-SDS012 


In  data  compression  classes  at  NMSU,  students  have  applied  various  compression 
techniques  to  data.  One  student  even  added  errors  to  see  the  effect  errors  would  have  on 
a  compressed  file.  In  an  effort  to  see  how  well  the  compression  could  work,  the  results 
for  an  image  is  summarized  below  (3): 


Compression 

Algorithm 

% 

Compression 

Channel  Error 
rate 

Errors 

made 

Comments 

None 

0% 

1x10'^  ' 

1 

Can’t  notice  error 

None 

0% 

1x10-^ 

58 

Slight  speckling 

None 

0% 

lxl0‘- 

5022 

Lots  of  speckles 

Huffman 

12% 

1x10'*’ 

1 

Shift  in  image 

Huffman 

12% 

Ixir' 

42 

Distorted,  but 
recognizable 

Huffrnan 

12% 

1x10'^ 

4562 

Decoder  failed 

Ziv-Lempel 

3% 

IxlO'" 

1 

Can’t  notice  error 

Ziv-Lempel 

3% 

1x10*^ 

58 

Highly  distorted 

Ziv-Lempel 

3% 

1x10'- 

4981 

Decoder  failed 

Arithmetic 

13% 

1x10’" 

1 

Decoder  failed 

Arithmetic 

13% 

“bdF 

58 

Decoder  failed 

Arithmetic 

13% 

IxTo^ 

4505 

Decoder  failed 

Rice 

22% 

1x10'" 

1 

Decoder  failed 

Rice 

22% 

1x10-^ 

5 

Decoder  failed 

Rice 

22% 

1x10'" 

32 

Decoder  failed 

From  this  table,  it  is  obvious  that  single  errors  with  compressed  data  can  be  disastrous.  It 
can  also  be  seen  that  when  the  Rice  algorithm  is  used  for  images  (for  which  it  was 
designed),  the  compression  is  almost  doubled.  If  the  Rice  algorithm  can  be  adapted  to 
handle  10,  12,  16  and  24  bits  per  word,  and  handle  variable  word  sizes  per  frame  (40, 

128,  160,  200,  256,  300  etc)  then  the  correlations  of  the  data  could  be  used  to  help  reduce 
the  file  size.  Since  the  data  from  one  frame  to  another  would  be  similar,  taking  advantage 
of  these  similarities  should  aid  in  the  compression  of  the  data.  One  way  the  Rice 
algorithm  uses  to  reduce  data  correlation  is  to  subtract  the  following  symbol  from  the 
previous  one.  This  works  well  for  images,  but  not  telemetry.  In  telemetry  an  entire 
frame  would  need  to  be  subtracted  from  the  following  frames  to  reduce  the  correlation. 
The  Rice  algorithm  is  not  currently  set  up  to  handle  the  variable  frame  sizes  needed. 

Channel  Error  Simulations 


Error  Sample  File  Description 

Error  files  were  obtained  from  Bert  Carner  at  Synthesys  Research  (4)  through  ARTM  for 
channel  error  simulations. 

General  Conditions: 

Frequency; L  band 
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Purpose  of  test:  A/B  comparison  of  FM  and  FQPSK-B  @  1  Mb/s 
TX  power:  nine  Watts,  each  transmitter 

TX  antenna  pattern:  well  behaved  omni  patterns  on  belly  of  aircraft,  single  antenna 
RX  antenna:  8-foot  dish,  bldg  5790 

Receivers:  Microdyne  model  700,  narrowband,  IF  BW  1  MHz,  single  channel,  i.e.,  NO 
COMBINING. 

Sample  point  conditions: 

Straight  and  level  flight,  Black  Mountain  corridor,  5k  foot  pressure  altitude, 
approximately  2500  foot  AGL,  200  Kt  indicated  airspeed 

The  two  points  represent  classic  multipath  events.  In  this  case  there  are  at  least  three  rays. 
There  is  a  strong  short  delay  ray  (10-20  nsec  range)  and  a  weak  long  delay  influence  (in 
the  range  of  2-5  usee) 

The  main  event  in  each  file  lasts  about  1.5-2  seconds  and  the  null  depth  is  deep,  at  least 
30  dB.  The  slant  range  is  only  about  40  miles,  very  strong  signal  conditions  (Eb/No 
greater  than  20  dB). 

Two  representative  files  were  used  to  demonstrate  the  effect  of  the  errors  on  the 
compressed  files.  For  the  data  files  SDS008,  and  SDS007  we  have  the  following  results: 


BERA  Simulation  Results 


File  name 

Compressed? 

Num  Errors 

Num  Bits 

BER 

Sds008.dat 

No 

3 

1105437696 

2.7138576-009 

Sds008.dat 

Yes 

2212390 

8650752 

0-2557 

Sds007.dat 

No 

10 

2.3475e+009 

4.25996-009 

Sds007.dat 

Yes 

218 

2.3475e+009 

9.28656-008 

The  compression  algorithm  used  for  compression  on  the  above  was  the  Ziv-Lempel 
algorithm  (since  it  tends  to  give  the  best  compression  results).  Here,  the  errors  can  cause 
major  problems.  Why  does  it  create  so  many  more  errors  in  the  compressed  case?  That 
is  based  on  how  the  individual  methods  of  compression  work.  With  Ziv-Lempel  (a 
dictionary  based  scheme),  an  error  in  the  data,  means  the  error  will  appear  in  the 
dictionary  and  so  affect  other  code  words.  So,  in  the  case  of  SDS008,  the  original  three 
errors  could  expand  to  2212390  errors  as  the  errors  go  through  the  dictionary  formed  by 
the  compression  algorithm.  However,  if  channel  coding  is  used,  these  problems  can  be 
avoided.  Other  ways  to  prevent  the  propagation  of  errors  is  to  reset  the  code  after  so 
many  bits  and  restart  the  code. 

Channel  Coding 

Looking  at  the  SDS008  case  where  3  errors  occur  in  the  original  file,  channel  coding 
would  be  able  to  eliminate  any  errors.  A  BCH  code  (1023,  973)  can  correct  5  errors  (5). 
The  973  is  the  number  of  incoming  bits  to  be  channel  coded  and  the  1023  is  the  number 
of  bits  sent.  For  each  973  bits  input  into  the  coder,  the  coder  will  add  50  bits  so  that  it 
can  find  up  to  and  including  5  errors  per  1023  bits  received.  For  this  file,  8891  blocks 
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will  need  to  be  channel  coded,  giving  a  total  of  9,095,302  bits  to  be  added  back  into  the 
file  to  eliminate  the  errors  caused  by  the  channel.  Still  at  this  rate,  the  amount  of 
compression  would  be: 

Compression  ratio  =  {1-  (8,650,752+9,095, 302)/(l, 105, 437,696) }  xlOO  =  98% 

This  still  gives  a  sizeable  amount  of  compression,  and  ensures  protection  from  the 
channel.  There  are  many  different  types  of  channel  noise  and  each  can  treated  using 
different  techniques.  BYU  is  currently  working  of  ways  of  handling  the  channel  effects 
that  are  pertinent  to  the  problems  encountered  by  the  military. 

Conclusions 


Data  compression  is  a  viable  option  to  aid  in  reducing  or  maintaining  bandwidth.  Even 
with  channel  errors  figured  in,  the  compression  achieved  is  worth  the  effort  involved. 
Significant  compression  can  be  achieved  through  off  the  shelf  technology.  Even  more 
compression  should  be  obtainable  with  modifications  of  the  Rice  algorithm. 
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Computer  Programs  Appendix 


James  Woods  wrote  a  lot  of  C  code  to  examine  the  statistics  of  the  data  files  and  lay  the 
preliminary  work  for  the  compression  algorithms.  The  data  and  the  initial  data  reading  programs 
were  obtained  from  Vernon  Diekmann,  at  TYBRIN  a  subcontractor  to  Edward’s  Air  Force  Base. 
Some  of  the  programs  written  by  Vem  Diekmann  were  modified  slightly  and  are  included  in  the 
list  of  programs  below.  The  compression  algorithms  used  were  fi'om  Mark  Nelson’s  book, 
except  for^tbe  Rice  algorithm,  which  was  obtained  through  NASA  and  the  University  of  New 
Mexico.  A  li.sting  of  the  programs  written  by  James  Woods  and  a  brief  description  of  what  they 
do  is  given  below: 

COMPRICE=>  This  program  allows  you  to  compare  a  file  that  has  been  run  through  rice  to 
another  file. 

COMP2  ,  COMP6  =>  This  program  just  compares  two  files  on  your  hard  drive. 

COMP3CD  =>  Takes  a  file  from  the  CD  and  compares  it  to  one  from  the  hard  drive. 

COUNT12  =>  Counts  12  words  configurations  from  the  hard  drive. 

FRAMAV  =>  Finds  the  average  frame,  the  most  popular  frame. 

FRAMINU1^>  Takes  the  first  frame  and  minuses  it  from  the  other  frames. 

FRAMMINU=>Takes  a  frame  that  was  saved  in  another  file  and  minuses  it  from  the  data  file. 
IMPROVE=>Same  thing  as  display  but  from  the  hard  drive. 

NEGATIVE=>Minuses  the  frames  without  negatives. 

NEW  TWO=>Counts  the  different  configurations  of  4  and  8  bits  from  the  hard  drive. 

READ  l-> Will  let  you  read  any  data  file  you  want  to  put  in. 

TEN  GREA=>Counts  the  different  configuration  of  10  bits  from  the  CD. 

TWOONE=>Counts  the  different  configurations  of  4  and  8  bit  from  the  hard  drive. 
TWOONECD=>Counts  the  different  configurations  of  4  and  8  from  the  CD. 

WORD=>Counts  the  different  configurations  of  12  bits  from  the  hard  drive. 

WORDCD=>Counts  the  different  configurations  of  12  bits  from  the  CD. 

The  FRAM  files  were  to  see  if  using  the  Rice  algorithm  where  the  frames  are  subtracted  would 
be  viable.  We  ran  into  some  limits  of  the  students’  programming  abilities  at  this  point  and  didn’t 
successfully  recreate  the  files  we  were  working  with  -  this  was  due  to  mapping  the  negative  data 
values  into  the  positive  values.  These  programs  should  not  be  run  without  fixing  the  mapping 
problem. 

The  Two  and  Word  programs  were  used  to  obtain  the  probability  densities  and  then  the  entropies 
of  the  individual  files. 

The  program  files  are  on  the  following  pages. 
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COMPRICE 


#include  <stdio.h> 

#include  <string.h> 

#define  numfiles  10 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles] [lengfiln],  ofilename [lengfiln],  ifilename[lengfiln],  outdat [lengfiln] 
plot[pltln+l]; 

char  fsll  [numfiles]  [lengfiln]; 

int  result,  status,  subframe,  fi,  a,  si^ei,  size2,  z; 

unsigned  int  data,  word,  interval,  rs,  sfid; 

unsigned  int  i,nwpf[numfiles],  nbpw [numfiles],  sfidw[numfiles],  sfido[numfiles], 

sff[numfiles],  sfl[numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  nfrms[numfiles],  wow[lengfiln]; 
unsigned  long  int  j==0,nf^0,  fsp [numfiles],  nbps [numfiles]; 
char  response,  excel[32],nev^le[lengfiln],  newfilel  [lengfiln]; 
char  hello; 

struct  tbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes] 

unsigned  int  nl  :  8; 

}; 

struct  kbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  ;  2; 
unsigned  int  b9  ;  2; 
unsigned  int  b  10:  2; 
unsigned  int  bl  1 :  2; 
unsigned  int  bl2;  2; 
unsigned  int  bl3:  2; 
unsigned  int  bl4:  2; 
unsigned  int  bl5:  2; 


unsigned  int  bl6:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  ;  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 


2048*/ 


struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8  = 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 


2048*/ 


union 


struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 


2048*/ 


union 


struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 


union 


struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 


struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 


2048*/ 


}  in3  ; 
union 

struct  Itbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible  [4096] ; 
struct  Ibytes  byte  [4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fin,  *fo,  *lp,  *conip;/*File  pointers*/ 
void  main() 


fin=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl  [i]  ,&fsll  [i]  ,&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 


&filename[i],&fsp[i],&nwpfIi],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl  [i]  ,&fsll  [i]  ,&nfrms[i])) 

{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 

fclose  (fin);/*Closed  cinput.dat*/ 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n'');/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpf[fi],nbpw[fi]); 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sffIfi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotai  Number  of  Frames  =  %ld\n", 
nbp  s  [fi]  ,nfrms  [fi] ) ; 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sfflfi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 
newfile[0]='e'; 
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newfile[l 

z=0; 

do 

{ 

newfiIe[2+z]=filename[fi][z];/*This  puts  'e:'  in  front  of  the  file 
name  so  it  can  read  it  off  the  cd*/ 

Z++; 

}while(z<=10);z=0; 

f=fopen(newfiIe,"rb");/*  Open  selected  data  file  */ 
printf("\nWhat  file  would  you  like  to  compare  that  to: "); 
scanf("%s",newfile  1 ); 

comp=fopen(newfilel,"rb");/*Open  compare  file  and  rb  mean  read  in  binary*/ 

printf("\nWhat  is  the  frame  size  after  rice:  "); 

scanf("%d",&sizel); 

printf("\nWhat  is  the  actual  frame  size  before  rice:  "); 
scanf("%d",&size2); 

printf("\n\n WORKING! ! !  !\n");/*Know  if  it  is  working*/ 

while  (ffead(in.data,(nwpf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1  ;/* Frame  counter*/ 


for(j=0;  j<nwpf[fi]*nbpw[fi]/8;  j-H-)/* Looping  through  the  frame*/ 

{ 

fread(inl.data.sizeof(int),l, comp);/* Read  in  data 

from  compared  file  and  store  it  in  ini. data*/ 

if(inl  .data[0]  !=in.data[j])/* Checking  if  the  two  files 

match  8  bits  at  a  time*/ 

{ 

if(in  1 .  data[  1  ] !  =in.  data|j  ] ) 

{ 

printf("\nERROR  at  frame  %ld  and 

in  %d\n",nf,j); 

printf("\nActual  %d  and  New  data 

%d\n",in.data[j],inl  .data[0]); 

printf("\nDo  you  want  to  return  to 

the  menu  (y  or  n):  "); 

response=getchar(); 

if(response=='y') 

{ 

goto  again; 

} 

} 

} 

} 

fread(in2.data.(size  1  -size2),  1  ,comp); 

} 

printf("\nDONE  COMPARING  %s  AND  %s\n\n",filename[fi],newfile); 
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printf  ("Successful  EndOfFile  reached  %Id  Minor  Frames.\n\n",  nf);/*  Cleanup,  close-up, 
and  start  over  */ 

fclose(f);/*Close  the  data  file*/ 
fclose(comp);/* Close  the  compare  file*/ 
nf  =  0; 

again;/*This  is  where  that  goto  function  goes  to*/ 

for  (i=0;  i  <  numfile;  i-i-i-){  printf("%d  %s\n",i,filename[i]);  }/*This  prints  menu  of  the 

files*/ 

printf  ("  Enter  the  file  Number\n"); 
printf  v"  or  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 


22 


C0MP2,  C0MP6 


include  <stdio.h> 

#include  <string.h> 

#define  numfiles  1 0 
#defme  lengfiln  40 
#define  pltln  125 

char  filename[numfiles][lengfiln],  ofilename[lengfiln],  ifilename [lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 

char  fsll[numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data,  word,  interval,  rs,  sfid; 

unsigned  int  i,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiies],  sfido[numfiles], 

sff[numfiles],  sfl[numfiies],  numfile,  fsl [numfiles],  nb; 
unsigned  int  nfrms [numfiles],  wow[lengfiln]; 
unsigned  long  int  j=0,nf=0,  fsp [numfiles],  nbps [numfiles]; 
char  response,  excel[32],newfile[numfiles] [lengfiln]; 
char  hello; 

struct  tbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles  { 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  b34es{ 

unsigned  int  nl  ;  8; 

}; 

struct  Itbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  ;  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  ;  2; 
unsigned  int  b9  ;  2; 
unsigned  int  b  10:  2; 
unsigned  int  bl  1:  2; 
unsigned  int  bl2:  2; 
unsigned  i^t  bl3:  2; 
unsigned  int  bl4:  2; 
unsigned  int  bl 5:  2; 


unsigned  int  bl6;  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in  ; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  ini  ; 
union 

struct  Itbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fin,  *fo,  *lp,  *comp;/*File  pointers*/ 

void  mainQ 

{ 

i  =  0; 

fm=fopen("cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 


«S:filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sficlw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i])) 

{  if  (i  <  numfiles-1)  i  =  i  +1; 

} 

numfile  =  i ; 

fclose  (fin);/*Closed  cinput.dat*/ 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  - 1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpf[fi],nbpw[fi]); 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],nfrms[fi]); 

printf("Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 

f=fopen(&filename[fi][0],"rb");/*  Open  selected  data  file  */ 
printf("\nWhat  file  would  you  like  to  compare  that  to:  "); 
scanf("%s",newfile[l]); 
hello=getchar();/*Take  up  the  return*/ 

comp=fopen(&newfile[l][0],"rb");/*Open  compare  file  and  rb  mean  read  in 

binary*/ 

response='f; 

printf("\nWORKING!!!!\n");/*Know  if  it  is  working*/ 

while  (fread(in.data,(nwpf];fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1  ;/*Frame  counter*/ 

if  (response  =  'f  ||  response  ==  'F')/*  Display  each  minor  frame  */ 

while  (fread(inl.data,(nwpf[fi]*nbpw[fi]/8),l,comp))/*Read  in 
data  from  compared  file  and  store  it  in  ini. data*/ 
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frame*/ 

match  8  bits  at  a  time*/ 
%ld",nfj); 


for(j=0;  j<nwpf[fi]*nbpw[fi]/8;  j-H-)/* Looping  through  the 

{ 

if(inl.data[j]!=in.data|j])/*Checking  if  the  two  files 

{ 

printf("\nERROR  in  frame  %ld  and  in 


printf("\n%02x\t%02x\n",in.dataO],inl .data[j]);/*Print  in  hex*/ 

printf("\nDo  you  want  to  return  to  menu 
(Y,N)\n’');/*If  the  is  a  error  this  will  let  you  return  to  the  begining*/ 

hello=getchar(); 

if(hello=='y'||hello=='Y') 

{ 

goto  again;/*  If  you  want  to  know 

where  this  goes  to  look  for  again  at  the  bottom*/ 

} 

} 

}break;/*This  breaks  the  while  loop  for  the  compare  file  so 

that  it  read  in  a  new  frame*/ 

} 

} 

printf(’'\nDONE  COMPARING  %s  AND  %s\n\n",filename[fi],&newfile[l][0]); 
printf  ("Successful  EndOfFile  reached  %ld  Minor  Frames.\n\n",  nf);/*  Cleanup,  close-up, 
and  start  over  */ 

fclose(f);/* Close  the  data  file*/ 
fclose(comp);/* Close  the  compare  file*/ 
nf  =  0; 

again:/*This  is  where  that  goto  function  goes  to*/ 

for  (i=0;  i  <  numfile;  i4-+){  printf("%d  %s\n",i,filename[i]);  }/*This  prints  menu  of  the 

files*/ 

printf  ("  Enter  the  file  Number\n"); 
printf  ("  or  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main*/ 


#include  <stdio.h> 

#include  <string.h> 

#defme  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles][lengfiln],  ofilename[lengfiln],  ifilename [lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 
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char  fsll  [numfiles]  [lengfiln] ; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data,  word,  interval,  rs,  sfid; 

unsigned  int  z=0,i,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido[numfiles], 
sff[numfiles],  sfl [numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  nfrms[numfiles],  wow[lengfiln]; 
unsigned  long  int  j=0,nf=0,  fsp[numfiles],  nbps[numfiles]; 
char  response,  excel[32],ne'wfile[lengfiln],file[lengfiln]; 
char  hello; 

struct  tbits 


{  - 

unsigned  int  bl  :  2; 
unsigned  int  b2  ;  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  :  8; 

}; 

struct  Itbits 

{ 

unsigned  int  bl  ;  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  ;  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  blO:  2; 
unsigned  int  bl  1:  2; 
unsigned  int  bl2:  2; 
unsigned  int  bl3:  2; 
unsigned  int  b  14:  2; 
unsigned  int  b  15:  2; 
unsigned  int  bl6:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 


unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  ;  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  ;  8; 

}; 

union 

struct  tbits  tbit[2048] ;/*  Max  allowable  Class  II  bytes  per  frame  1 63  84/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  b>^e[2048]; 
unsigned  char  data[2048]; 

}  ini  ; 
union 

struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte[4096]; 
unsigned  long  int  da1a[4096]; 

}  out ; 

FILE  *f,  *fm,  *fo,  *lp,  *comp:/*File  pointers*/ 

void  mainQ 

{ 

i  =  0; 

fm=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl  [i]  ,&fsll  [i]  ,&nfrms  [i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpfIi],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i])) 

{  if  (i  <  numfiles-1)  i  =  i  +1; 
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printf("\nDo  you  want  to  return  to  menu 

(Y  N)\n");/*If  the  is  a  error  this  will  let  you  return  to  the  begining*/ 

hello=getchar();/*Take  up  the  return*/ 
if(hello=’y'||hello=='Y’) 

{ 

goto  again;/*  If  you  want  to  know 

where  this  goes  to  look  for  again  at  the  bottom*/ 

}break;/*This  breaks  the  while  loop  for  the  compare  file  so 

that  it  read  in  a  new  frame*/ 

} 

printf("\nDONE  COMPARING  %s  AND  %s\n\n",filename[fi],file);/*This  prints  the  file 
names  that  were  compared*/ 

printf  ("Successful  EndOfFile  reached  %ld  Minor  Frames. \n\n",  nf);/*  Cleanup,  close-up, 

and  start  over  */ 

fclose(f);/* Close  the  data  file*/ 
fclose(comp);/*Close  the  compare  file*/ 
nf=0; 

again:/*This  is  where  that  goto  function  goes  to*/ 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]);  }/*This  prints  menu  of  the 

files*/ 

printf  ("  Enter  the  file  Number\n"); 
printf  ("  or  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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C0MP3CD 


#include  <stdio.h> 

#include  <string.h> 

#define  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles] [lengfiln],  ofilename[lengfiln],  ifilename [lengfiln],  outdat[lengfiln], 
plot[pltIn+l]; 

char  fsll  [numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data,  word,  interval,  rs,  sfid; 

unsigned  int  z=0,i,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido[numfiles], 

sff[numfiles],  sflfnumfiles],  numfile,  fsl[numfiles],  nb; 
unsigned  int  nfrms [numfiles],  wow[lengfiln]; 
unsigned  long  int  j=0,nf=0,  fsp[numfiles],  nbps[numfiles]; 
char  response,  excel[32],newfile[lengfiln],file[lengfiln]; 
char  hello; 

struct  tbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  :  8; 

}; 

struct  kbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  b  10:  2; 
unsigned  int  bl  1:  2; 
unsigned  int  bl2:  2; 
unsigned  int  bl3:  2; 
unsigned  int  b  14:  2; 
unsigned  int  bl5:  2; 


unsigned  int  bl6:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in  ; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte  [2048]; 
unsigned  char  data[2048]; 

}  ini  ; 
union 

struct  Itbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 
unsigned  long  int  data[4096]; 

}  out; 

FILE  *f,  *fm,  *fo,  *lp,  *comp;/*File  pointers*/ 

void  mainQ 

{ 

i  =  0; 

fm=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 


&filename[i],&fsp[i],&nwpfti],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i])) 

{  if  (i  <  numfiIes-1)  i  =  i +1; 

} 

numfile  =  i ; 

fclose  (fin);/*Closed  cinput.dat*/ 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpf[fi],nbpw[fi]); 

printf("Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi]  ,sfido[fi]  ,sfftfi]  ,sfl  [fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %Id\n", 
nbps  [fi]  ,nfrms[fi] ) ; 

printf("Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sff[f])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 
newfile[0]-e'; 
newfile[l 
z=0; 
do 
{ 

newfile[2+z]=filename[fi][z];/*This  puts  'e:'  in  front  of  the  file 
name  so  it  can  read  it  off  the  cd*/ 

Z++; 

}  while(z<=  1 0);z=0 ; 

f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 

printf("\nType  in  the  path  and  file  name  you  would  like  to  compare  that  to:\n"); 
scanf("%s",file); 

hello=getchar();/*Take  up  the  return*/ 

comp=fopen(file,"rb");/*Open  compare  file  and  rb  mean  read  in  binary*/ 
printf("\n\nWORKING! ! !  !\n\n");/*Know  if  it  is  working*/ 


while  (fread(in.data,(nwpf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1  ;/*  Frame  counter*/ 

while  (fread(inl.data,(nwpfIfi]*nbpw[fi]/8),l,comp))/*Read  in 
data  from  compared  file  and  store  it  in  ini. data*/ 

for(j=0;  j<nwpf[fi]*nbpw[fi]/8;  j++)/*Looping  through  the 


frame*/ 

match  8  bits  at  a  time*/ 


{ 


if(inl.dataIj]!=in.data[j])/*Checking  if  the  two  files 

%ld",nf,j);/*Print  error  if  the  data  doesn't  match*/ 

printf("\n%02x\t%02x\n",in.dataO],inl  .dataO]);/* Print  in  hex*/ 


printf("\nERROR  in  frame  %ld  and  in 


printf("\nDo  you  want  to  return  to  menu 

(Y,N)\n");/*If  the  is  a  error  this  will  let  you  return  to  the  begining*/ 

hello=getchar();/*Take  up  the  return*/ 
if(hello=’y'l|hello=='Y') 

{ 

goto  again;/*  If  you  want  to  know 

where  this  goes  to  look  for  again  at  the  bottom*/ 

} 

} 

}break;/*This  breaks  the  while  loop  for  the  compare  file  so 

that  it  read  in  a  new  frame*/ 

} 

printf("\nDONE  COMPARING  %a  AND  %s\n\n’',filename[fi],file);/*This  prints  the  file 
names  that  were  compared*/ 

printf  ("Successful  EndOfFile  reached  %ld  Miner  rrames.\n\n",  nf);/*  Cleanup,  close-up, 
and  start  over  */ 

fclose(f);/*  Close  the  data  file*/ 
fclose(comp);/*Close  the  compare  file*/ 
nf=0; 

again:/*This  is  where  that  goto  function  goes  to*/ 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]);  }/*This  prints  menu  of  the 

files*/ 

printf  ("  Enter  the  file  Number\n"); 
printf  ("  or  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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C0UNT12 


#include  <stdio.h> 

#include  <string.h> 

#defme  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles] [lengfiln],  ofilename[Iengfiln],  ifilename[lengfiln],  outdat[lengfiln] 
plot[pItln+l]; 

char  fsll  [numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data, word,  interval,  rs,  sfid; 

unsigned  int  b=0,i,  j,  z=0,x=0,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles], 
sfido[numfiles], 

sff[numfiles],  sfl [numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  n=0,  nfhns [numfiles],  wow[lengfiln]; 
unsigned  long  int  nf=0,fsp [numfiles], 

nbps'[numfiles],afbl[256],counterl[256],afb2[256],counter2[256],afb3[256],counter3[256] 
char  response,  excel[32]; 

char  hello,newfilel[30],newfile2[30],newfile3[30],  newfile[40]; 
struct  tbits 
{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

\- 
) » 

struct  bytes  { 

unsigned  int  nl  :  8; 

}; 

struct  Itbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  intb9  : 2; 
unsigned  int  blO:  2; 
unsigned  int  bll:  2; 
unsigned  int  b  12:  2; 
unsigned  int  bl3:  2; 


unsigned  int  bl4;  2; 
unsigned  int  bl5:  2; 
unsigned  int  bl6:  2; 

}; 

struct  Iniblesj 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

{ 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  b>tes  per  frame  1 6384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in  ; 
union 
{ 

struct  Itbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  -  4096*/ 

"  struct  Inibles  nible[4096]; 

struct  Ibytes  byte[4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fin,  *fo,  *lpl,  *lp2,  *lp3; 

void  main() 

{i  =  0; 

fm=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw'[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl  [i]  ,&fsll  [i]  ,&nfrms  [i] )) 


36 


{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fin); 


for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1 ;  result=scanf("%d’',&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi]  ,fsp[fi]  ,nwpf[fi]  ,nbpw[fi] ); 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],(Iong)  nfrms[fi]); 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb-0;i=l; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 

printfC'What  is  the  file  you  want  counted:\n"); 
scanf("%s",newfile); 
printf("\nWORKING! ! !  !\n"); 


zero*/ 


zero*/ 


f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 
n-O; 
z-0; 
do 
{ 

afb2[z]=z;/*This  puts  numbers  into  the  array  for  comparing*/ 
counter2[z]=0;/*This  puts  zeros  in  the  array  to  set  the  count  to 

Z++; 

}while(z<=15);z=0;/*Puts  numbers  up  to  15  in  the  array*/ 
do 
{ 

afbl  [n]=n;/*This  puts  numbers  into  the  array  for  comparing*/ 
counterl[n]=0;/*This  puts  zeros  in  the  array  to  set  the  count  to 


n++; 

}while(n<=255);n=0;/*Puts  numbers  up  to  255  in  the  array*/ 
while  (fread(in.data,(nwpf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1  ;/* Frame  counter*/ 

j=0; 
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match*/ 


match*/ 


bit  data  into  nibles*/ 


frame  of  12  bit  words*/ 


for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j++) 

{ 

do 

{ 

if(in.data[j]==afbl[n])/*Check  if  the  numbers 

{ 

counter  l[n]+=l;/* counter  if  they 
break;/* Breaks  once  counted*/ 

} 

n++; 

}while(n<=255);n=0; 

} 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  12 

{  out.nible[i]  .n3  =  in.nible|j]  .n2; 
out.nible[i]  .n2  =  in.nible|j]  .nl; 
out.nible[i]  .nl  =  in.nible[j+l].n2; 
out.nible[i+l].n3  =  in.nible[j+l].nl; 
out.nible[i+l].n2  =  in.nible[j+2].n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 
i=i+2; 

} 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/12;  j++)/*Loop  through  the 
{  z=0; 

do 


{ 

if(out.nible[j].nl==afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nible[j].n2=afb2[z])/*  Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nible[j].n3==afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 


} 

Z++; 

}while(z<=15);z=0; 
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} 

printf("\n%s  DONE ! ! !  !\n",newfile); 
printf("\nOutput  file  name  for  the  two  hex  #: "); 
scanf("%s",newfilel); 

printf("\nOutput  file  name  for  the  one  hex  #: "); 
scanf("%s",newfile2); 

lp2=fopen(newfile2,''w");/*Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lp2,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
lpl=fopen(newfilel,''w");/*Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lpl,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
n=0; 
do 

fprintf(lpl,"%lx\tCount\t%ld\n",afbl [n],(Iong)counterl  [n]);/* Printing  to  the 

output  file*/ 

n++; 

}while(n<=255);n=0; 

z=0; 

do 

fprintf(lp2,"%lx\tCount\t%ld\n",afb2[z],(long)counter2[z]);/*Printing  to  the 

output  file*/ 

Z++; 

}while(z<=l  5);z=0; 

printf  ("\nSuccessful  EndOfFile  reached  %ld  Minor  Frames.\n\n",(long)  nf);/*  Cleanup, 
close-up,  and  start  over  */ 

fclose(f); 
fclose(lpl); 
fclose(lp2); 
nf  =  0; 

for  (i=0;  i  <  numfile;  i+-i-){  printf("%d  %s\n",i,filename[i]);  }/*?rints  menu  of  files*/ 
printf  ("  Enter  the  file  Number  ;  \n"); 
printf  ("  <ctrl>  c  to  end;  \n"); 

} 

}  /*  End  Main*/ 
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FRAMAV 


#include  <stdio.h> 

#include  <string.h> 

#define  numfiles  1 0 
#deflne  lengfiln  40 
#define  pltln  125 

char  filename[numfiles][lengfiln],  ofilename[lengfiln],  ifilename[lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 

char  fsll  [numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data, word,  interval,  rs,  sfid; 

unsigned  int  z=0,  i,  j,c=0,  nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido[numfiles], 
sff[numfiles],  sfl[numfiles],  numfile,  fsl[numfiles],  nb; 
unsigned  int  wow[lengfiln]; 

unsigned  long  int  nfrms[numfiles],nf=0,  k=0,  fsp [numfiles],  nbps[numfiles],  abfilengfiln], 

abfl^lengfiln],  abf2 [lengfiln]; 

char  response,  excel[32],  outfile[40],newfile[40]; 

char  hello; 

struct  tbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  ;  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles  { 

unsigned  int  nl  ;  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  :  8; 

}; 

struct  kbits 


unsigned  int  bl  ;  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  ;  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  b  10:  2; 
unsigned  int  bll:  2; 
unsigned  int  bl2:  2; 
unsigned  int  bl3:  2; 
unsigned  int  bl4:  2; 


unsigned  int  bl5;  2; 
unsigned  int  bl6:  2; 


}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  ;  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

^  struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in ; 
union 

^  struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8  = 

2048*/ 

struct  nibles  mble[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  ini  ; 
union 

^  struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}in2; 

union 

^  struct  tbits  tbit[2048];/*  Max  allowable  Class  II  b\'tes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte  [2048]; 


unsigned  char  data[2048]; 


}  in3  ; 
union 

struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fin,  *fo,  *lp,  '^fbut; 

void  main() 

{ 

for  (i  =  0;  i  <  pltln;  i++)  plot[i]  = '  ';/*  Initialize  array  */ 
i  =  0; 

fin=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename  [i]  ,&fsp  [i]  ,&nwpf[i]  ,&nbpw[i]  ,&nbps  [i]  ,&sfidw[i]  ,«&sfido  [i]  ,&sff[i]  ,&sfl  [i]  ,& 
fsl[i],&fsll[i],&nfnns[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl  [i]  ,&fsll  [i],&nfnns[i])) 

{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fin); 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  - 1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharQ;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename  [fi]  ,fsp[fi]  ,nwpf[fi]  ,nbpw[fi]); 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %Id\tTotal  Number  of  Frames  =%ld 
\n",nbps[fi],nfrms[fi]); 

printf("Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 
fsl[fi],fsll[fi]);  ' 
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nb=0;i=l; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 
newfile[0]='e'; 
newfile[l]=';'; 
z=0;  * 

do 

{ 

newfile[2+z]=filename[fi][z]; 

Z++; 

}  while(z<=  1 0)  ;z=0 ; 
printf("\nWORKING! ! !  !\n"); 
f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 
for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j++) 

{ 

abflj]=0; 

} 

while  (fread(in.data,(nwpf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1 ; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j++) 

{ 

abf|j]=abf|j]+in.dataO]; 

} 

} 

printfC'Type  in  path  and  file  name  for  the  output:\n"); 
scanf("%s",outfile); 
fout=fopen(outfile,"w"); 
for  (k=0;  k<nwpf[fi]*nbpw[fi]/8;  k++) 

{ 

abfl[k]=abf[k]/nf; 
fprintf(fout,"%ld",abfl  [k]); 

} 

printf  ("Successful  EndOfFile  reached  %d  Minor  Frames.\n",  nf);/*  Cleanup,  close-up, 
and  start  over  */ 
fclose(f); 
fclose(fout); 
if  (result  >  2)fclose(fo); 
nf  =  0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i, filename [i]); } 
printf  ("  Enter  the  file  Number  :  \n"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main*/ 
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FRAMINUl 


#include  <stdio.h> 

#include  <math.h> 

#include  <string.h> 

#define  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles] [lengfiln],  ofilename[lengfiln],  ifilename [lengfiln],  outdat[lengfiln], 
pIot[pltln+l]; 

char  fsll[nunifiles][lengfiln]; 

int  result,  status,  subframe,  cool,  abfl  [lengfiln],  abf4[lengfiln],  abfZ [lengfiln],  bye,  a,  u; 
unsigned  int  data, word,  interval,  rs,  sfid,  joke; 

unsigned  int  good,  k=0,i=0,  z=0,  j,  nwp^numfiles],  nbpw[numfiles],  sfidw[numfiles], 
sfido[numfiles], 

sff[numfiles],  sfl [numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  b=0,  wow[lengfiln],  abf3 [lengfiln]; 

unsigned  long  int  nfrms[numfiles],  nf=0,  abf5[500],  fsp [numfiles],  nbps [numfiles]; 
long'^ int  boo=0,  one,  negative[300],  negcountpOO],  h,  m,  p,  positivepOO],  poscount[300]; 
char  response,  check,  okay,  excel[32],  newfile[40],  file[40],  addfile[40],  could[40]; 
char  hello,  song; 

long  int  1=0,  count=0,  c=0,  v,  saw,  hi,  fi,  try,  abf[lengfiln]; 
struct  tbits 
{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes] 

unsigned  int  nl  :  8; 

}; 

struct  kbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  blO:  2; 
unsigned  int  bll:  2; 


unsigned  int  bl2:  2; 
unsigned  int  bl3:  2; 
unsigned  int  b  14:  2; 
unsigned  int  bl5:  2; 
unsigned  int  bl6:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

struct  tbits  tbit[6048] ;/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
unsigned  char  data[6048]; 

}  in  ; 
union 

struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
signed  char  data[6048]; 

}in4; 

union 

struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
unsigned  char  data[6048]; 

}  ini  ; 
union 
{ 


2048*/ 


struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 


struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
unsigned  char  data[6048]; 

}in2; 

union 

struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

■2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
signed  char  data[6048]; 

}  in3  ; 
union 

struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
signed  char  data[6048]; 

}  in5  ; 
union 
{ 

struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  -  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *g,  *fm,  *fo,  *lp,  *fout,  *added,  *frame; 

void  main() 

{ 

printf("\nMENU:"); 

printf("\nADD->'a"'); 

printf("\nMINUS->'m"'); 

printf("\nAVG->’v’\n:"); 

response=getchar() ; 

hello=getchar(); 

check=response; 

okay=response; 

i  =  0; 
c=0; 

count=0; 

1=0; 

fin=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fm,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 


&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sffIi],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  —  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n  , 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sffIi],&sfl[i],& 

fsl  [i]  ,&fsll  [i]  ,&nfrms[i])) 

{  if  (i  <  numflles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fin); 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf(’'%d",&fi);  fi  !=  - 1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi]  ,fsp  [fi]  ,nwpflfi]  ,nbpw[fi] ); 

printfC’Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sffIfi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],nfrms[fi]); 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=I; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 

1=0; 

if(response=='m'l|response=-'M') 

{ 

newfile[0]='e’; 
newfile[l  ]=':'; 
b=0; 
do 

newfile[2+b]=filename[fi][b];/*This  puts  'e:'  in  front  of  the  file 

name  so  it  can  read  it  off  the  cd*/ 

b++; 

}  while(b<=  1 0)  ;b=0; 

f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 

fo=fopen("frame.dat","wb"); 

fout=fopen("  new.  dat" , "  wb " ) ; 

fread(inl  .data,(nwpf[fi]*nbpw[fi]),l  ,f); 

fwrite(inl.data,(nwpf[fi]*nbpw[fi]),l,fo); 
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fclose(fo); 

v=0; 

1=0; 

c=0; 

count=l; 

while(fread(&in2 .  data[0] ,  sizeof(int),  1  ,f)) 

{ 

in4.data[0]=(inl.data[c]-in2.data[0])*(-l); 

fwrite(&in4.data[0],sizeof(in4.data[0]Xl^fout); 

iix.data[0]=inl.data[c]+in4.data[0]; 

if(in4.data[0]<0) 

{ 

do 

{ 

if(negati  ve  [m]  ==in4 .  data  [0] ) 

{ 

negcount[m]+=l; 

h++; 

} 

I++; 

}  while(l<=m);l=0; 
if(h!=l) 

{ 

negative  [m+ 1  ] =in4 .  data[0] ; 
negcount[m]+=l; 

} 

h=0; 

} 

if(in4.data[0]>=0) 

{ 

do 

{ 

if(positive[p]  ==in4 .  data[0]) 

{ 

poscount[p]+=l; 

h++; 

} 

1++; 

}  while(l<=p);l=0; 
if(h!=l) 

{ 

positi  ve[p+ 1  ] =in4.data[0] ; 
poscount[p]+=l; 

} 

h=0; 

} 

if(in.data[0]  !=in2.data[0]) 

{ 
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printf("(Error)  in  %d  =  ini  %d  +  in4 

%d\n",in.data[0],inl.data[c],in4.data[0]); 

hello=getchar(); 

} 

if(count=(nwpf[fi]  *nbpw[fi]/8)) 

{ 

count=l; 

1++; 

} 

C++; 

count++; 

if(c==(nwpflfi]  *  nbpw[fi])) 

{ 

c=0; 

printfC’."); 

} 

in5.data[0]=(inl.data[c]-in2.data[l  ])*(-!); 
fwrite(&in5.data[0],sizeof(in5.data[0]),l,fout); 
in.data[0]=inl.data[c]+in5.data[0]; 
if(in.data[0]  !=in2.data[l  ]) 

{ 

printf("(Error)  in  %d  =  ini  %d  +  in4 

%d\n",in.data[0],inl.data[c],in4.data[0]); 

hello=getchar(); 

} 

if(count==(nwpf[fi]*nbpw[fi]/8)) 

{ 

count=l; 

1++; 

} 

C++; 

count++; 

if(c==(nwpf[fi]*nbpw[fi])) 

{ 

c=0; 

printfC'."); 

} 

hello=getchar(); 

} 

printf("\nThe  minus  output  files  are  frame.dat  and  new.dat  in  this 

directory  \n\n"); 

printf  ("Successful  EndOfFile  reached  %ld  Minor  Frames.\n",(long)l); 

fclose(f); 

fclose(fout); 

count=0; 

} 

1=0; 

c=0; 

boo=0; 
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if(response=='a'|  Iresponse— A') 

{ 

fclose(f); 

newfile[0]='e'; 

newfile[l]=':'; 

b=0; 

do 

{ 

newfile[2+b]=filename[fi][b];/*This  puts  'e:'  in  front  of  the  file 
name  so  it  can  read  it  off  the  cd*/ 

b++; 

}while(b<=10);b=0; 
f^fopen(ne  wfile,  "rb " ) ; 

printf("\nWhat  is  the  file  you  what  to  add  to:"); 
scanf("%s", could); 

printf("\nOutput  file  name  for  the  added  data: "); 

scanf("%s",addfile); 

hello=getchar(); 

g=fopen(could,"rb ") ; 

added=fopen(addfile,"wb"); 

fi-ame==fopen("frame.dat","rb"); 

fread(in.data,(nwpf[fi]*nbpw[fi]),l , frame); 

printf("\nWas  the  frames  minused:  (y,n) "); 

song=getchar(); 

hello=getchar(); 

if(song=='y') 

{ 

fwrite(in.data,(nwpf[fi]*nbpw[fi]),  1  .added); 
fread(in5.data,(nwp^fi]*nbpvv[fi]),  1  ,f); 

} 

fclose(frame); 

1=0; 

c=l; 

while(fread(&in3  .data[0]  ,sizeof(in3  .data[0] ),  1  ,g)) 

{ 

in2.data[0]=in3.data[0]+in.data[l]; 
fread(&inl  .data[0],sizeof(inl  .data[0]),l,f); 
if(in2.data[0]  !=inl  .data[0]) 

{ 

printf("frame  %ld\n",c); 
printf("word  %d\n",l); 
printf("\n(Error)  ini  %d  in2 

%d\n",inl  .data[0],in2.data[0]); 

printf("\nIN3  %d  +  in  %d  =  in2 
%d\n",in3.data[0],in.data[l],in2.data[0]); 

hello=getchar(); 

} 

1++; 

if(l==(nwpf[fi]  *nbpw[fi])) 
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{ 

printfC’."); 

1=0; 

C-H-; 

} 

fwrite(&in2.data[0],sizeof(in2.byte[0]),l, added); 

} 

printf("\nThe  added  output  file  %s  is  done!\n",addfile); 

printf("\nNumber  of  frames  %ld\n",(long)c); 

fclose(added); 

fclose(g); 

fclose(f); 

fclose(frame); 

} 

1=0; 

c=0; 

if(response==’v'||response=='v') 

{ 

newfile[0]-e'; 

newfile[l]- 

file[0]='e'; 

file[l]=':'; 

b=0; 

do 

newfile[2+b]=filename[fi][b];/*This  puts  'e;'  in  front  of  the  file 

name  so  it  can  read  it  off  the  cd*/ 

file[2+b]=filename[fi][b]; 

b++; 

}while(b<=10);b=0; 

f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 

fo=fopen("frame.dat","wb"); 

fout=fopen("new.dat","wb"); 

k=0; 

j=0; 

while(fread(in  1  .data,sizeof(in.byte[0]),  1  ,f)) 

{ 

ifG<=l) 

{ 

abf[k]=inl  .data[0]; 

} 

ifa>i) 

{ 

abf[k]=inl  .data[0]+abf[k]; 

} 

k-H-; 

if(k=(nwpfl^fi]  *  nbpw[fi])) 

{ 

printfC'."); 
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} 


} 


k=0; 

j++; 


do 

{ 

in2.data[v]=abf[v]/j ; 

V-H-; 

}  while(v<=(nwpflfi]  *nbpw[fi])) ;  v=0 ; 
fvvrite(in2.data,(nwpf[fi]  *nbpw[fi]),  1  ,fo); 
c=0; 
k=0; 

fclose(fo); 

fclose(f); 

g=fopen(file,"rb"); 

while(fread(&in.data[0]  ,sizeof(int),  1  ,g)) 

{ 

in4.data[0]=(in2.data[k]-in.data[0])*(-l); 
fwrite(&in4.data[0],sizeof(in4.data[0]),  1  ,fout); 
k++; 

in5.data[0]=(in2.data[k]-in.data[l])*(-l); 
fwrite(&in5.data[0],sizeof(in5.data[0]),l  .fout); 
k-H-; 

if(k==(nwpf|^fi]  *nbpw[fi]/8)) 

{ 

1++; 

} 

if(k==(nwpf[fi]  *nbpw[fi])) 

{  printfC’."); 
k=0; 

} 

} 

fclose(g); 

fclose(fout); 

fclose(frame); 

fclose(f); 

fclose(fo); 

printf("\nThe  average  frame  in  is  frame.dat  and  the  minus  data  is  in 

new.dat  \n\n"); 

printf  ("\nSuccessful  EndOfFile  reached  %ld  Minor 
Frames.\n",(long)l);/*  Cleanup,  close-up,  and  start  over  */ 

} 

else 

{} 

k=0; 

boo=0; 

j-0; 

c=0; 
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nf  =  0; 

1=0; 

printf("\nMENU:"); 

printf("\nADD->’a'"); 

printf("\nMINUS->'m"'); 

printf("\nAVG->’v'\n:"); 

scanf("%c'',&response); 

hello=getchar(); 

count=0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  ros'\n",i,fiIename[i]);  }/*Print  out  menu*/ 
printf  ("  Enter  the  file  Number  :  \n"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main*/ 
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FRAMMINU 


#include  <stdio.h> 

#include  <math.h> 

#include  <string.h> 

#define  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename  [numfiles]  [lengfiln],  ofilename[lengfiln],  ifilename[lengfiln],  outdat[lengfiln] 
plot[pltln+l]; 

char  fsll[numfiles]  [lengfiln]; 

int  result,  status,  subframe,  cool,  abfl  [lengfiln],  abf4 [lengfiln],  abl2 [lengfiln],  bye,  a,  u; 
unsigned  int  data, word,  interval,  rs,  sfid,  joke; 

unsigned  int  good,  k=0,i=0,  z=0,  j,  nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles], 
sfido[numfiles], 

sff[numfiles],  sfl [numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  b=0,  wow[lengfiln],  abO [lengfiln]; 

unsigned  long  int  nfrms [numfiles],  nf=0,  abf5[1000],  fsp [numfiles],  nbps[numfiles]; 
long  int  boo=0,  one; 

char  response,  check,  okay,  excel[32],  newflle[40],  file[40],  addfile[40],  could[40]; 
char  hello,  song; 

long  int  1=0,  count=0,  c=0,  v,  saw,  hi,  fi,  try,  abf[lengfiln]; 
struct  tbits 
{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  ;  2; 
unsigned  int  b4  ;  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  ;  8; 

}; 

struct  kbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  ;  2; 
unsigned  int  b9  :  2; 
unsigned  int  blO:  2; 
unsigned  int  bl  1:  2; 


unsigned  int  bl2;  2; 
unsigned  int  bl3:  2; 
unsigned  int  bl4:  2; 
unsigned  int  bl5:  2; 
unsigned  int  b  16:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

^  struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
unsigned  char  data[6048]; 

}  in  ; 
union 

^  struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
signed  char  data[6048]; 

}  in4; 
union 

struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte  [6048]; 
unsigned  char  data[6048]; 

}  ini  ; 
union 
{ 


struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8  = 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte[6048]; 
unsigned  char  data[6048]; 

}  in2  ; 
union 

struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte  [6048]; 
signed  char  data[6048]; 

}in3; 

union 

struct  tbits  tbit[6048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[6048]; 
struct  bytes  byte  [6048]; 
signed  char  data[6048]; 

}  in5  ;  ^ 

union 

{ 

struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte[4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *g,  *fm,  *fo,  *lp,  *fout,  *  added,  *  frame; 

void  main() 

{ 

printf("\nMENU;'’); 

printf("\nADD->’a'"); 

printf("\nMINUS->'m"'); 

printf("\nAVG->V\n:"); 

response=getchar(); 

hello=getchar(); 

check=response; 

okay=response; 

i  =  0; 
c=0; 

count=0; 

1=0; 

fm=fopen("e:cinput.dat",'’rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fm,''%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 
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&filename  [i]  ,&fsp  [i]  ,&nwpf{i],&nbpw[i]  ,&nbps[i]  ,&sfidw[i]  ,&sfido[i]  ,&sff[i]  ,&sfl  [i]  ,& 
fsl[i]  ,&fsll[i]  ,&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpfIi],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sffIi],&sf][i],& 

fsl  [i]  ,&fsll  [i]  ,&nfrms[i])) 

{  if  (i  <  numfiles-1)  i  =  i  +1; 

} 

numfile  =  i ; 
fclose  (fm); 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpf[fi],nbpvv[fi]); 

printfC’Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %Id\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],nfrms[fi]); 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=I; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l ; 

1=0; 

if(response=='m'llresponse=-'M') 

{ 

newfile[0]=’e'; 

newfile[l]=';'; 

b=0; 

do 

{ 

newfile[2+b]=filename[fi][b];/*This  puts  'e:'  in  front  of  the  file 
name  so  it  can  read  it  off  the  cd*/ 

b++; 

}while(b<=10);b=0; 

f=fopen(newfile,"rb'');/*  Open  selected  data  file  */ 

fo=fopen( "  frame.dat" ,  "wb ") ; 
foub=fopen( "  new.dat" , "  wb " ) ; 
fread(inLdata,(nwpfi;fi]  *nbpw[fi]),  1  ,f); 
fwrite(inl  .data,(nwpflfi]*nbpw[fi]),  1  ,fo); 
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fclose(fo); 

v=0; 

1=0; 

c=0; 

count=l; 

while(fread(&in2.data[0],sizeof(int),  1  ,f)) 

in4.data[0]=(inl.data[c]-in2.data[0])*(-l); 
fwrite(&in4.data[0],sizeof(in4.data[0]),  1  ,fout); 
iii.data[0]=inl.data[c]+in4.data[0]; 
if(iii';data[0]  !=in2.data[0]) 

printf("(Error)  in  %d  =  ini  %d  +  in4 

%d\n",in.data[0],inl.data[c],in4.data[0]); 

helIo=getchar(); 

} 

if(count==(nwpf[fi]  *nbpw[fi]/8)) 

{ 

count=l; 

1++; 

} 

C++; 

count++; 

if(c==(nwpf[fi]  *nbpw[fi])) 

{ 

c=0; 

printfC'."); 

} 

in5  .data[0]=(in  1  .data[c]-in2.data[  1  ])*  (- 1 ); 
fwrite(&in5  .data[0]  ,sizeof(in5  .data[0] ),  1  ,fout); 
iii.data[0]=in  1  .data[c]+in5  .data[0] ; 
if(in.data[0]!=in2.data[l]) 

{ 

printf("(Error)  in  %d  =  ini  %d  +  in4 

%d\n",in.data[0],inl.data[c]  jn4.data[0]); 

hello=getchar(); 

} 

if(count=(nwpf[fi]  *nbpw[fi]/8)) 

{ 

count=l; 

1++; 

} 

C++; 

count++; 

if(c==(nwpf[fi]*nbpw[fi])) 

{ 

c=0; 

printfC'."); 

} 
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} 

printf("\nThe  minus  output  files  are  frame.dat  and  new.dat  in  this 

directory  \n\n"); 

printf  ("Successful  EndOfFile  reached  %ld  Minor  Frames.\n",(long)l); 

fclose(f); 

fclose(fout); 

count=0; 

} 

1=0; 

c=0; 

boo=0; 

if(response=='a'||response='A’) 

{ 

fclose(f); 

newfile[0]='e'; 

newfile[l]- 

b=0; 

do 

{ 

newfile[2+b]=filename[fi][b];/*This  puts  'e:'  in  front  of  the  file 
name  so  it  can  read  it  off  the  cd*/ 

b++; 

}while(b<=10);b=0; 

f=fopen(newfile,"rb"); 

printf("\nWhat  is  the  file  you  what  to  add  to:"); 
scanf("%s", could); 

printf("\nOutput  file  name  for  the  added  data: "); 

scanf("%s",addfile); 

hello=getchar(); 

g=fopen(could, "  rb " ) ; 

added=fopen(addfile, "  wb") ; 

frame=fopen("  frame .  dat" ,  "rb ") ; 

fread(in.  data, (nwpf[fi]  *nbpw[fi]),  1 , frame); 

printf("\nWas  the  frames  minused:  (y,n) "); 

song=getchar(); 

hello=getchar(); 

if(song==y) 

{ 

fwrite(in.data,(nwpflfi]*nbpw[fi]),l, added); 
fi-ead(in5.data.(nw'p^fi]*nbpw[fi]),l,f); 

} 

fclose(frarae); 

1=0; 

c=l; 

while(fread(&in3  .data[0],sizeof(in3  .data[0]),  1  ,g)) 

{ 

in2.data[0]=in3  .data[0]+in.data[l] ; 
fread(&inl  .data[0],sizeof(inl  .data[0]),  1  ,f); 
if(in2.data[0]  !=inl  .data[0]) 


{ 

printf("frame  %ld\n",c); 
printf("word  %d\n",l); 
printf("\n(EiTor)  ini  %d  in2 

%d\n",inl  .data[0],in2.data[0]); 

printf("\nIN3  %d  +  in  %d  =  in2 
%d\n",in3.data[0],in.data[l],in2.data[0]); 

hello=getchar(); 

} 

1++; 

if(l==(nwpf[fi]  *nbpw[fi])) 

{ 

printfC'."); 

1=0; 

C++; 

} 

fwrite(&in2.data[0],sizeof(in2.byte[0]),l, added); 

} 

printf("\nThe  added  output  file  %s  is  done!\n",addfile); 

printf("\nNumber  of  frames  %ld\n",(long)c); 

fclose(added); 

fclose(g); 

fclose(f); 

fclose(frame); 

} 

1=0; 

c=0; 

if(response=='v’||response=='v') 

{ 

newfile[0]='e'; 

newfile[l]=':'; 

file[0]='e’; 

file[l]=’:'; 

b=0; 

do 

{ 

newfile[2+b]=filename[fi][b];/*This  puts  'e:'  in  front  of  the  file 
name  so  it  can  read  it  off  the  cd*/ 

file[2+b]=filename[fi][b]; 

b++; 

}while(b<=l  0);b=0; 

f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 

fo=fopen("frame.dafV'wb"); 

fout=fopen( "  new.dat " , "  wb " ) ; 

k=0; 

j=0; 

while(fread(in  1  .data,sizeof(in.byte[0]),  1  ,f)) 

{ 

if(j<-l) 
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{ 


abf[k]=inl  .data[0]; 


} 


} 

ifG>i) 

{ 

abf[k]=in  1  .data[0]+abf[k] ; 

} 

k-H-; 

if(k==(nwpflfi]  *nbpw[fi])) 

{ 

printfC'."); 

k=0; 

j++; 


do 

{ 

in2.data[v]=abf[v]/j ; 


V-H-; 

}while(v<=(nwpf[fi]*nbpw[fi]));v=0; 
fwrite(in2.data,(nwpf[fi]  *nbpw[fi]),  1  ,fo); 
c=0; 
k=0; 

fclose(fo); 

fclose(f); 

g=fopen(file,"rb"); 

while(fread(&in.data[0],sizeof(int),  1  ,g)) 

{ 

in4.data[0]=(in2.data[k]-in.data[0])*  (- 1 ); 

fwrite(&in4.data[0],sizeof(in4.data[0]),l,fout); 

k++; 

in5 .  data[0]=(in2 .  data[k] -in.  data[  1  ] )  *  (- 1 ) ; 

fwrite(&in5.data[0],sizeof(in5.data[0]),l,fout); 

k++; 

if(k==(nwpf[fi]  *nbpw[fi]/8)) 

{ 

1++; 

} 

if(k==(nwpf[fi]  *  nbpw[fi] )) 

{  printfC'."); 
k=0; 

} 

} 

fclose(g); 

fclose(fout); 

fclose(frame); 

fclose(f); 

fclose(fo); 


printf("\nThe  average  frame  in  is  frame.dat  and  the  minus  data  is  in 


new.dat  \n\n"); 


printf  ("\nSuccessfuI  EndOfFile  reached  %ld  Minor 
Frames.\n",(long)l);/*  Cleanup,  close-up,  and  start  over  */ 


} 

else 

{} 

k=0; 

boo=0; 

j=0; 

c=0; 

nf=0: 

1-0; 

printf("\nMENU:"); 

printf("\nADD->’a"'); 

printf("\nMINUS->'m"’); 

printf("\nAVG->'v'\n:"); 

scanf("  %c " ,  &response) ; 

hello=getchar(); 

count=0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]);  }/*Print  out  menu*/ 
printf  ("  Enter  the  file  Number  :  \n"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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IMPROVE 


/* 

*- 


♦ 

♦ 

Name: 

Display  Datax 

* 

Creation  Date: 

09-Sep-97 

* 

Description: 

Display  Data  from  Telemetry  Sample  Data  Sets 

♦ 

Author: 

Vernon  Diekinann,  TYBRIN 

* 

* 

Target  Machine/OS/Compiler: 

Iiitel  Pentium/Windows  95/NT/MS  Visual  C-h- 

* 

* 

* 

Revision  History 

* 

* 

Initials  Date 

Description 

* 

VD 

09/09/97 

Initial  version 

* 

VD 

10/10/97 

Added  12  and  16  bit 

* 

VD 

01/15/98 

Added  24  bit  and  read  byte  records 

* 

VD 

01/28/98 

Input  the  file  name  for  the  data.dat  file. 

* 


*/ 

#include  <stdio.h> 

#include  <string.h> 

#defme  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles] [lengfiln],  ofilename[lengfiln],  ifilename[lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 

char  fsll[numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data, word,  interval,  rs,  sfid; 

unsigned  int  i,  j,  nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido [numfiles], 
sff[numfiles],  sfl[numfiles],  numfile,  fsl  [numfiles],  nb; 
unsigned  int  nfrms[numfiles],  wow[lengfiln]; 
unsigned  long  int  nf=0,  fsp [numfiles],  nbps [numfiles]; 
char  response,  excel[32]; 
char  hello; 

struct  tbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
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}; 

struct  bytes  { 

unsigned  int  nl  :  8; 

}; 

struct  kbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  b  10:  2; 
unsigned  int  bl  1:  2; 
unsigned  int  b  12:  2; 
unsigned  int  bl3:  2; 
unsigned  int  bl4:  2; 
unsigned  int  bl 5:  2; 
unsigned  int  b  16:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

{ 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in ; 
union 
{ 


struct  Itbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [40 96]; 
unsigned  long  int  data[4096]; 

}  out; 

FILE  *f,  *fin,  *fo,  *Ip; 

void  main() 

{ 

for  (i  =  0;  i  <  pltln;  i-H-)  plot[i]  = '  ';/*  Initialize  array  */ 
i  =  0; 

fin=fopen("cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename  [i]  ,&fsp[i]  ,&nwpf[i]  ,&nbpw[i]  ,&nbps[i]  ,&sfidw[i]  ,&sfido  [i]  ,&sff[i]  ,&sfl  [i]  ,& 
fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fm,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i]  ,&nwpf[i]  ,&nbpvv[i]  ,&nbps[i]  ,&sfidw[i]  ,&sfido[i]  ,&sff[i],&sfl  [i]  ,& 
fsl  [i]  ,&fsll  [i]  ,&nfrms[i] )) 

{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fin); 

for  (i=0;  i  <  numfile;  i++){  printf("%d-— >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf(''%d",&fi);  fi  !=  -1;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpf[fi],nbpw[fi]); 

printf("Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],nfrms[fi]); 

printf("Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l ; 

printf  ("  Enter\nF  to  print  each  frame  or  \nW  to  plot  a  single  word:  '\n"); 
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response  =  getchar(); 


if  (response  =  'w'  ||  response  ==  'W') 

{ 

printf  ("  Enter  the\n  Word,  Subframe,  Interval  to  display;  \n"); 
result  =  scanf("%d%d%d",&word,&subframe,&interval); 
hello  =  getchar(); 

printf  ("  Enter  'c'  to  continue  or  type  the  path  and  file  name  for  the  data  of 
the  word  selected:  \n"); 

result  =  scanf("%s",&excel); 

} 

f=fopen(&filename[fi][0],"rb");/*  Open  selected  datafile  */ 
result  =  strlen(&excel[0]); 
if(result  >  2) 

fo=fopen(&excel[0],"w");/*  Open  option  data  output  file  for  excel  */ 

while  (ffead(in.data,(nwpf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 

{  nf  =  nf  +  1 ; 

if  (response  =  'f  ||  response  ==  'F')/*  Display  each  minor  frame  */ 

{  printf  ("\nFrame=%d\n",nf); 

for  (j=0:  j<nwpf[fi]*nbpw[fi]/8;  j-H-) 

{ 

printf("%02x'',in.data|j]); 

} 

printf("\n"); 

printf  ("Enter  to  continue  <ctrl>C  to  quit"); 
hello=getchar(); 

} 

else  if  (response  ==  'w'  ||  response  ==  'W')  /*  Display  selected  data 

parameter  */ 

{ 

if  (nbpw[fi]==8) 

{ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=l)/*  Unravel  packed  8 

bit  data  */ 

{ 

out.byte[i].nl  =  in.bytelj].nl; 

} 

} 

if  (nbpw[fi]==10) 

{ 

i=0; 

■  for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=5)/*  Unravel  packed  10 

bit  data  */ 

{ 
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out.tbit[i].b5  =  in.tbit[j].b4; 
out.tbit[i].b4  =  in.tbit[j].b3; 
out.tbit[i].b3  =  in.tbit[j].b2; 
out.tbit[i].b2  =  in.tbit[j].bl; 
out.tbit[i].bl  =  in.tbit|j+l].b4; 
out.tbit[i+l].b5  =  in.tbit[j+l].b3; 
out.tbit[i+l].b4  =  in.tbit[j+l].b2; 
out.tbit[i+l].b3  =  in.tbitO+l]-bl; 
out.tbit[i+l].b2  =  in.tbit[j+2].b4; 
out.tbit[i+l].bl  =  in.tbit|j+2].b3; 
out.tbit[i+2].b5  =  in.tbitO+2].b2; 
out.tbit[i+2].b4  =  in.tbit[j+2].bl; 
out.tbit[i+2].b3  =  in.tbit[j+3].b4; 
out.tbit[i+2].b2  =  in.tbit[j+3].b3; 
out.tbit[i+2].bl  =  in.tbit[j+3].b2; 
out.tbit[i+3].b5  =  in.tbit[j+3].bl; 
out.tbit[i+3].b4  =  in.tbit[j+4].b4; 
out.tbit[i+3].b3  =  in.tbit[j+4].b3; 
out.tbit[i+3].b2  =  in.tbit[j+4].b2; 
out.tbit[i+3].bl  =  in.tbit[j+4].bl; 
i=i+4; 

} 

} 

if  (nbpw[fi]-=12) 

{ 

i=0; 

for  (i=0;  j<nwpf[fi]*nbp\v[fi]/8;  j+=3)/*  Unravel  packed  12 

bit  data  */ 

{  out.nible[i]  .n3  =  in.nible[j]  .n2; 
out.nible[i]  .n2  =  in.nible[j]  .nl; 
out.nible[i]  .nl  =  in.nible[j+l].n2; 
out.nible[i+l].n3  =  in.nible[j+l].nl ; 
out.nible[i+l].n2  =  in.nible[j+2].n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 
i=i+2; 

} 

} 

if  (nbpw[fi]==16) 

{ 

i-0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=2)/*  Unravel  packed  16 

bit  data  */ 

{  out.byte[i].n2  =  in.byte[j]  .nl; 
out.byte[i].nl  =  in.byte[j+l].nl; 
i=i+l; 

} 

} 

if  (nbpw[fi]==24) 
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bit  data  */ 


sub  frame  data  */ 


data*pltln; 


out.data[word-l]); 


commutated  data  */ 


data*pltln; 


out .  data[word- 1  +j  ] ) ; 


{ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  24 

{  out.byte[i].n3  =  in.bytelj]  .nl; 

out.byte[i].n2  =  in.byte[j+l].nl; 
out.byte[i].nl  =  in.byte[j+2].nl; 

H+1; 

} 

} 

if  ((sfidw[fi]  >=•  0)  &  (  interval  !=  0)  &  (subframe  >=  0))/*  Handle 

{  sfid  =  out.data[sfidw[fi]-l]; 
sfid  =  sfid  »  rs; 

for  (j=0;  j<(sfl[fi]-sff[fi]+l);  j+=interval) 

{  if  ((subframe+j)  ==  sfid) 

{  data  =  (unsigned)out.data[word-l]/pltln; 
data  =  (unsigned)out.data[word-l]- 

plot[data]  = 
plot[data+l]  =  '\0'; 
printf  ("%d%s\n",nf,  plot); 
plot[data]  = ' 
plot[data+l]  = ' 

if  (result  >  2)fprintf  (fo,"%d\t%d\n",nf, 

} 

} 

} 

else  if  ((subframe  =--!)&(  interval  !=  0))/*  Handle  super 

{  for  (j~0;  j<=nwpf[fi] ;  j+=interval) 

{  data  =  (unsigned)out.data[word-l+j]/pltln; 
data  =  (unsigned)out.data[word-l+j]- 

plot[data]  = 
plot[data+l]  =  '\0'; 
printf  ("%d%s\n",nf,  plot); 
plot[data]  = ' '; 
plot[data+l]  = ' '; 

if  (result  >  2  )fprintf(fo,"%d\t%d\n",nf, 

} 

} 

else  /*  Just  plain  old  data  */ 

{  data  =  (unsigned)out.data[word- 1  ]/pltln; 

data  =  (unsigned)out.data[word-l]-data*pltln; 
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plot[data]  = 
plot[data+l]  =  '\0'; 
printf  ("%d%s\n",nf,  plot); 
plot[data]  = 
plot[data+l]  = ' 

if  (result  >  2)fprintf  (fo,"%d\t%d\n",nf.  out.data[word-l]); 

} 

} 

else 

{  /*  Do  nothing  */ 

} 

} 

printf  ("Successful  EndOfFile  reached  %d  Minor  Frames.\n",  nf);/*  Cleanup,  close-up, 
and  start  over  */ 
fclose(f); 

if  (result  >  2)fclose(fo); 
nf  =  0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]); } 
printf  ("  Enter  the  file  Number  : 'm"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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NEGATIVE 


#include  <stdio.h> 

#include  <string.h> 

#defme  numfiles  10 
-#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles] [lengfiln],  ofilename[lengfiln],  ifilename[lengfiln],  outdat[lengfiln], 
plot[pItln+l]; 

char  fsll  [numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data, word,  interval,  rs,  sfid; 

unsigned  int  b=0,i,  j,  z=0,x=0,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles], 
sfido  [numfiles], 

sfl^numfiles],  sfl [numfiles],  numfile,  fsl[numfiles],  nb; 
unsigned  int  n=0,  nfrms[numfiles],  wow[lengfiln]; 
unsigned  long  int  nf=0,fsp [numfiles], 

nbps[numfiles],afbl[256],counterl[256],afb2[256],counter2[256],afb3[256],counter3[256]; 
char  response,  excel[32]; 

char  hello,newfilel[30],newfile2[30],newfile3[30],  newfile[40]; 
long  negativel  [2000],negative2[2000],negative3[2000],  negative4[2000]; 
struct  tbits 
{ 

signed  int  bl  :  2; 
signed  int  b2  :  2; 
signed  int  b3  :  2; 
signed  int  b4  :  2; 

}; 

struct  nibles{ 

signed  int  nl  :  4; 
signed  int  n2  :  4; 

}; 

struct  bytes  { 

signed  int  nl  :  8; 

}; 

struct  Itbits 

{ 

signed  int  bl  :  2; 
signed  int  b2  :  2; 
signed  int  b3  :  2; 
signed  int  b4  :  2; 
signed  int  b5  :  2; 
signed  int  b6  :  2; 
signed  int  b7  :  2; 
signed  int  b8  :  2; 
signed  int  b9  :  2; 
signed  int  blO:  2; 
signed  int  bl  1;  2; 
signed  int  b  12:  2; 


signed  int  bl3:  2; 
signed  int  bl4:  2; 
signed  int  b  15:  2; 
signed  int  b  16:  2; 

}; 

struct  lnibles{ 

signed  int  nl  :  4; 
signed  int  n2  :  4; 
signed  int  n3  :  4; 
signed  int  n4  :  4; 
signed  int  n5  :  4; 
signed  int  n6  :  4; 
signed  int  n7  :  4; 
signed  int  n8  :  4; 

}; 

struct  lbytes{ 

signed  int  nl  :  8; 
signed  int  n2  :  8; 
signed  int  n3  :  8; 
signed  int  n4  :  8; 

}; 

union 

{ 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048], 
struct  bytes  byte[2048]; 
signed  char  data[2048]; 

}  in ; 
union 
{ 

struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 
signed  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fin,  *fo,  *lpl,  *lp2,  *Ip3; 

void  mainQ 
{i  =  0; 

fin=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fm,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sffIi],&sfl[i],& 
fsl  [i]  ,&fsll  [i]  ,&nfims  [i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %Id  %d  %d  %d  %d  %d  %s  %ld\n", 
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&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i]  [i]  ,&nfrms  [i])) 

{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  -  i ; 
fclose  (fin); 

for  (i^O;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  G;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1 ;  result=scanf("%d'',&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharQ;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi]  ,fsp  [fi]  ,nwpf[fi]  ,nbpw[fi] ) ; 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi]  ,sfido[fi]  ,sff[fi]  ,sfl  [fi]) ; 

printffNumber  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],(long)  nfrms[fi]); 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 


zero*/ 


zero*/ 


printf("What  is  the  file  you  want  countedAn"); 
scanf("%s",newfile); 
printr("\nWORKING! ! !  !\n"); 

f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 
n=0; 
z=0; 
do 

^  afb2[z]=z;/*This  puts  numbers  into  the  array  for  comparing*/ 
counter2[z]=0;/*This  puts  zeros  in  the  array  to  set  the  count  to 


Z++; 

}while(z<=15);z=0;/*Puts  numbers  up  to  15  in  the  array*/ 
do 

afbl[n]=n;/*This  puts  numbers  into  the  array  for  comparing*/ 
counterl  [n]=0;/*This  puts  zeros  in  the  array  to  set  the  count  to 


n++; 

}while(n<=255);n=0;/*Puts  numbers  up  to  255  in  the  array*/ 
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match*/ 

match*/ 


while  (fread(in.data,(nwpf{fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1  ;/*  Frame  counter*/ 

j=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j++) 

{ 

do 

if(in.dataO]==afbl  [n])/*Check  if  the  numbers 

I 

counterl  [n]+=l  ;/*counter  if  they 
break;/*Breaks  once  counted*/ 

} 


}while(n<=255);n-0; 

} 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j++) 

{ 

do 

{ 

if(in.data|j]==(-l  *afbl  [n])) 

{ 

negative  l[n]+=l; 
printf("%l  count 

%r',in.data|j],negativel  [n]); 

hello=getchar(); 

break; 

} 

n++; 

}whilt(n<=255);n=0; 

} 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  12' 

bit  data  into  nibles*/ 

{  out.nible[i]  .n3  =  in.nible[j]  .n2; 
out.nible[i]  .n2  =  in.nible[j]  .nl; 
out.nible[i]  .nl  =  in.nible[j+l].n2; 
out.nible[i+l].n3  =  in.nible[j+l].nl; 
out.nible[i+l].n2  =  in.nible[j+2].n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 

i=i+2; 

} 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/12;  j++)/*Loop  through  the 

frame  of  12  bit  words*/ 

{z=0; 

do 

{ 
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if(out.nible|j].nl=afb2[z])/*  Checks  if  the 
niimbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nible[j].n2=afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nibIe[i].n3^  afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 


frame  of  12  bit  words*/ 


%1"  ,out.  nible  [j  ]  .n  1  ,negative2  [n] ) ; 


%l",out.nible|j].n2,negative2[n]); 


%l",out.nible[j].n2,negative2[n]); 


Z++; 

}while(z<=15);z=0; 

} 

z=0; 

for  (i=0;  j<nwpf[fi]*nbpw[fi]/12;  j++)/*Loop  through  the 

{  z=0; 

do 

{ 

if(out.nible[j].nl==(afb2[z]*(-l))) 

{ 

negative2[z]+=l; 
printf("%l  count 

hello=getcbar(); 

} 

if(out.  nible  [j  ]  .n2==(afb2  [z]  *  (- 1 ))) 

{ 

negative2[z]+=l; 
printf("%l  count 

heIlo=getchar(); 

} 

if(out.nible(j  ]  .n3==(afb2  [z]  *  (- 1 ))) 

{ 

negative2[z]+=l; 
printf("%l  count 

hello=getchar(); 

} 

Z++; 

}  while(z<=l  5);z=0; 

}  printfC'."); 
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} 

printf("\n%s  DONE! ! ! !\n",newfile); 
printf("\nOutput  file  name  for  the  two  hex  #:  "); 
scanf("%s",newfilel); 

printf("\nOutput  file  name  for  the  one  hex  #: "); 
scanf("%s",newfile2); 

lp2=fopen(newfile2,"w");/* Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lp2,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 

Ip l=fopen(newfilel,"w");/* Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lpl,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
n=0; 
do 

fprintf(lpl,"%lx\tCount\t%ld\n",afbl[n],(long)counterl[n]);/*Printing  to  the 

output  file*/ 

n++; 

}while(n<=255);n=0; 

z=0; 

do 

fprintf(lpl,"%lx\tCount\t%ld\n",(afbl  [n]  *(-!)), (long)negativel  [n]); 
n-H-; 

}while(n<=255);n=0; 

z=0; 

do 

fprintf(lp2,"%lx\tCount\t%ld\n",afb2[z],(long)counter2[z]);/*Printing  to  the 

output  file*/ 

Z++; 

}while(z<=15);z=0; 

do 

{ 

fprintf(lp2,"%lx\tCount\t%ld\n",(afbl[n]*(-l)),(long)negative2[n]); 

Z++; 

}while(z<=15);z=0; 

printf  ("\nSuccessful  EndOfFile  reached  %ld  Minor  Frames.\n\n",(long)  nf);/*  Cleanup, 
close-up,  and  start  over  */ 

fclose(f); 

fcloseOpl); 

fclose(lp2); 
nf  =  0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]);  }/*Prints  menu  of  files*/ 
printf  ("  Enter  the  file  Number  :  \n''); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 


75 


76 


NEW  TWO 


#include  <stdio.h> 

#include  <string.h> 

#define  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles] [lengfiln],  ofilenanie[lengfiIn],  ifilename[lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 

char  fsll  [numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data, word,  interval,  rs,  sfid; 

unsigned  int  b=0,i,  J,  z=0,x=0,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles], 
sfido[numfiles], 

sff[numfiles],  sfl [numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  n=0,  nfrms[numfiles],  wow[lengfiln]; 
unsigned  long  int  nf=0,fsp [numfiles], 

nbps  [numfiles]  ,afb  1  [25 7]  ,counter  1  [257]  ,afb2  [25  7]  ,counter2  [257]  ,afb3  [2  5  7]  ,counter3  [257]; 
char  response,  excel[32]; 

char  hello,newfilel[30],newfile2[30],newfile3[30],  newfile[40]; 
struct  tbits 
{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  ;  2; 
unsigned  int  b4  :  2; 

}; 

struct  niblesf 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  :  8; 

}; 

struct  kbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  ;  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  blO:  2; 
unsigned  intbll:  2; 
unsigned  int  b  12:  2; 
unsigned  int  bl3:  2; 


unsigned  int  bl4:  2; 
unsigned  int  bl5:  2; 
unsigned  int  bl6:  2; 


}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4;  f  * 

unsigned  int  n7  ;  4; 
unsigned  int  n8  ;  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in  ; 
union 

struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  “  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 

unsigned  long  int  data[4096]; 

}  out; 

unsigned  long  int  temp; 

FILE  *f,  *fm,  *fo,  *lpl,  *lp2,  *lp3; 

void  main() 

{i  =  0; 

fm=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fm,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sffti],<S:sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 
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&fllename  [i]  ,&fsp[i]  ,&nwpf[i]  ,&nbpw[i]  ,&nbps[i]  ,&sfidw[i]  ,&sfido[i]  ,&sff[i]  ,&sfl  [i]  ,& 
fsl[i],&fsll[i],&nfrms[i])) 

{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fin); 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  ;  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpf[fi],nbpw[fi]); 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],(long)  nfrms[fi]); 

printf("Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 

newfile[0]='e'; 

newfile[l]=':'; 

b=0; 

do 

{ 

newfile[2+b]=filename[fi][b];/*This  puts  'e:'  in  front  of  the  file 
name  so  it  can  read  it  off  a  cd*/ 

b++; 

}  while(b<=  1 0)  ;b=0 ; 
printf("\nWORKING! ! !  !\n"); 


f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 
n=0; 
z=0; 
do 

{ 

afb2[z]=z;/*This  puts  numbers  into  the  array  for  comparing*/ 
counter2[z]=0;/*This  puts  zeros  in  the  array  to  set  the  count  to 

zero*/ 

Z-H-; 

}while(z<=15);z=0;/*Puts  numbers  up  to  15  in  the  array*/ 
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zero*/ 


afbl[n]=n;/*This  puts  numbers  into  the  array  for  comparing*/ 
counterl  [n]=0;/*This  puts  zeros  in  the  array  to  set  the  count  to 

n-H-; 

}while(n<=255);n=0;/*Puts  numbers  up  to  255  in  the  array*/ 
while  (fread(in.data,(nv/pf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1  ;/* Frame  counter*/ 

j=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j-H-) 

{ 

do 


match*/ 


match*/ 


bit  data  into  nibles*/ 


frame  of  12  bit  words*/ 


{ 

if(in.dataO]==afbl[n])/*Check  if  the  numbers 

{ 

counterl  [n]+=l  ;/*counter  if  they 
break;/*Breaks  once  counted*/ 

} 

n++; 

}while(n<=255);n=0; 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  12 

{  out.nible[i]  .n3  =  in.nible[i]  .n2; 
out.nible[i]  .n2  =  in.niblelj]  .nl; 
out.nible[i]  .nl  =  in.nible[j+l].n2; 
out.nible[i+l].n3  =  in.nible[j+l].nl; 
out.nible[i+l].n2  =  in.nible[j+2].n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 
i=i+2; 

} 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/12;  j++)/*Loop  through  the 
{z=0; 

do 


{ 

if(out.nible[i].nl=afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nible[i].n2==afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l ; 
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} 

if(out.nible[j].n3=afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

Z“Hh; 

}while(z<=15);z=0; 


} 

printf("\n%s  DONE! ! !  !\n",filename[fi]); 
printf("\nOutput  file  name  for  the  two  hex  #: "); 
scanf("%s",newfile  1 ); 

printf("\nOutput  file  name  for  the  one  hex  #: "); 
scanf("%s "  ,newfile2); 

lp2=fopen(newfile2,"w");/*Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lp2,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
lpl=fopen(newfilel,’'w");/*Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lpl,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
n=0; 
do 

{ 

if(counterl  [n]<counterl  [n-1]) 

{  printfC'..."); 

temp=counter  1  [n- 1  ] ; 

counter  1  [n-l]=counterl  [n]; 

counterl[n]=temp; 

afbl[257]-afbl[n-l]; 

afbl[n-l]=afbl[n]; 

afb  1  [n]=afb  1  [2  5  7] ; 

} 

n++; 

}while(n<=255);n=0; 

do 

fprintf(lpl,"%lx\tCount\t%ld\n",afbl[n],(long)counterl[n]);/*Printing  to  the 

output  file*/ 

n++; 

}while(n<=255);n=0; 

z=0; 

do 

^  fprintf(lp2,"%lx\tCount\t%ld\n",afb2[z],(long)counter2[z]);/*Printing  to  the 

output  file*/ 

Z++; 

}  while(z<=l  5);z=0; 

printf  ("\nSuccessful  EndOfFile  reached  %ld  Minor  Frames.\n\n",(long)  nf);/*  Cleanup, 
close-up,  and  start  over  */ 
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fclose(f); 
fclose(lpl); 
fclose(lp2); 
nf  =  0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]);  }/*Prints  menu  of  files*/ 
printf  ("  Enter  the  file  Number  ;  \n"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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READl 


/* 

*. 

* 

* 

* 

* 

♦ 

* 

*. 

* 

* 


Name:  Display_Data.c 

Creation  Date:  09-Sep-97 

Description:  Display  Data  from  Telemetry  Sample  Data  Sets 

Author:  Vernon  Diekmann,  TYBRIN 

Target  Machine/OS/Compiler:  Intel  Pentium/Windows  95/NT/MS  Visual  C++  5.0 


Revision  History 


*  Initials  Date  Description 

*  _  _ 


* 

VD 

09/09/97 

* 

VD 

10/10/97 

* 

VD 

01/15/98 

* 

VD 

01/28/98 

Initial  version 

Added  12  and  16  bit 

Added  24  bit  and  read  byte  records 

Input  the  file  name  for  the  data.dat  file 


*/ 

#include  <stdio.h> 

#include  <string.h> 

#defme  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles][lengfiln],  ofilename[lengfiln],  ifilename[lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 

char  fsll[numfiles][lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data, word,  interval,  rs,  sfid; 

unsigned  int  i,  j,  nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido [numfiles], 
sff[numfiles],  sfl [numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  nf=0,  nfrms [numfiles],  wow[lengfiln]; 
unsigned  long  int  fsp [numfiles],  nbps [numfiles]; 
char  response,  excel  [3  2]: 
char  hello; 

struct  tbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
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}; 

struct  bytes  { 

unsigned  int  nl  ;  8; 

}; 

struct  Itbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

unsigned  int  b5  :  2;  , 

unsigned  int  b6  :  2; 

unsigned  int  b7  :  2; 

unsigned  int  b8  :  2; 

unsigned  int  b9  :  2; 

unsigned  int  blO;  2; 

unsigned  int  bll:  2; 

unsigned  int  b  12:  2; 

unsigned  int  bl3:  2; 

unsigned  int  bl4:  2; 

unsigned  int  bl5:  2; 

unsigned  int  bl6:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

^  struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in; 

union 

{ 
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struct  Itbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte[4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fm,  *fo,  *lp; 

void  main() 

{ 

•c  X  ■ 

for  (i  =  OiJ  <  pltln;  i++)  plot[i]  =  ";/*  Initialize  array  */ 


i  =  0; 

fin=fopen("cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fm,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %d\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfiri],& 
fsl  [i]  ,&fsll  [i]  ,&nfrms[i])) 

{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fin); 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d’',&fi);  fi  !=  -1;  result-scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpf[fi],nbpw[fi]); 

printfC’Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],nfrms[fi]); 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fl]-nb+l; 

printf  ("  Enter\nF  to  print  each  frame  or^nW  to  plot  a  single  word:  \n"); 
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response  =  getcharQ; 


if  (response  =  'w'  ||  response  =  'W') 

{ 

printf  ("  Enter  the\n  Word,  Subframe,  Interval  to  display:  \n"); 
result  =  scanf("%d%d%d",&word,&subframe,&interval); 
hello  =  getcharO; 

printf  ("  Enter  'c'  to  continue  or  type  the  path  and  file  name  for  the  data  of 
the  word  selected:  \n"); 

result  =  scanf("%s",&excel); 

} 

f=fopen(filename[fi],"rb");/*  Open  selected  data  file  */ 
result  =  strlen(&excel[0]); 
if(result  >  2) 

fo=fopen(&excel[0],"w");/*  Open  option  data  output  file  for  excel  */ 

while  (fread(in.data,(nwpf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 

{  nf  =  nf  +  1; 

if  (response  ==  'f  ||  response  ==  'F')/*  Display  each  minor  frame  */ 

{  printf  ("\nFrame=%d\n",nf); 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j++) 

{ 

printf("%02x",in.data[j]); 

} 

printf("\n"); 

printf  ("Enter  to  continue  <ctrl>C  to  quit"); 
hello=getchar(); 

} 

else  if  (response  ==  'w'  ||  response  ==  'W')  /*  Display  selected  data 

parameter  */ 

{ 

if  (nbpw[fi]==8) 

{ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=l)/*  Unravel  packed  8 

bit  data  */ 

{ 

out.byte[i].nl  =  in.byte|j].nl; 

} 

} 

if  (nbpw[fi]==10) 

{ 

i-0; 

for  (i=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=5)/*  Unravel  packed  10 

bit  data  */ 

{ 
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bit  data  */ 


bit  data  */ 


out.tbit[i].b5  = 
out.tbit[i].b4  = 
out.tbit[i].b3  = 
out.tbit[i].b2  = 
out.tbit[i].bl  = 
out.tbit[i+l].b5 
out.tbit[i+l].b4 
out.tbit[i+l].b3 
out.tbit[i+l].b2 
out.tbit[i+l].bl 
out.tbit[i+2].b5 
out.tbit[i+2].b4 
out.tbit[i+2].b3 
out.tbit[i+2].b2 
out.tbit[i+2].bl 
out.tbit[i+3].b5 
out.tbit[i+3].b4 
out.tbit[i+3].b3 
out.tbit[i+3].b2 
out.tbit[i+3].bl 
i=i+4; 

} 

} 


in.tbit|j].b4; 

in.tbit[j].b3; 

in.tbit[j].b2; 

in.tbit[j].bl; 

in.tbit[j+l].b4; 

=  in.tbitO+l]-b3; 
=  in.tbit[j+l].b2; 
=  in.tbit[j+l].bl; 
=  in.tbit|j+2].b4; 
=  in.tbit|j+2].b3; 
=  in.tbitO+2].b2; 
=  in.tbit[j+2].bl; 
=  in.tbit[j+3].b4; 
=  in.tbit[j+3].b3; 
=  in.tbit|j+3].b2; 
=  in.tbit[j+3].bl; 
=  in.tbitO+4].b4; 
=  in.tbit[j+4].b3; 
=  in.tbit[j+4].b2; 
=  in.tbit[j+4].bl; 


if  (nbpw[fi]==12) 

{ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  12 

{  out.nible[i]  .n3  =  in.nible[j]  .n2; 
out.nible[i]  .n2  =  in.nible[j]  .nl; 
out.nible[i]  .nl  =  in.nible|j+l].n2; 
out.nible[i+l].n3  =  in.nible[i+l].nl; 
out.nible[i+l].n2  =  in.nible[j+2].n2; 
out. nible[i+ lj.nl  =  in.nible[j+2].nl; 
i=i+2; 

} 

} 

if  (nbpw[fi]==16) 

{ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=2)/*  Unravel  packed  16 

{  out.byte[i].n2  =  in.byteO]  .nl; 
out.byte[i].nl  =  in.byte[j+l].nl; 
i=i+ 1 ; 

} 

} 

if  (nbpw[fi]— 24) 
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{ 


bit  data  */ 


sub  frame  data  */ 


data*pltln; 


out.data[word-l]); 


commutated  data  */ 


data*pltln; 


out.data[word-l+j]); 


i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  24 

{  out.byte[i].n3  =  in.byteQ]  .nl; 

out.byte[i].n2  =  in.byteQ+lJ.nl; 
out.b5de[i].nl  =  in.byte[j+2].nl; 
i=i+l; 

} 

} 

if  ((sfidw[fi]  >=  0)  &  ( interval  !=  0)  &  (subframe  >=  0))/*  Handle 

{  sfid  =  out.data[sfidw[fi]-l]; 
sfid  =  sfid  »  rs; 

for  (j=0;  j<(sfl[fi]-sffi;fi]+l);  j+=interval) 

{  if  ((subframe+j)  ==  sfid) 

{  data  =  (unsigned)out.data[word- 1  ]/pltln; 
data  =  (unsigned)out.data[word-l]- 

plot[data] 
plot[data+l]  =  ’\0’; 
printf  ("%d%s\n",nf,  plot); 
plot[data]  =  "; 
plot[data+l]  = ' 

if  (result  >  2)fprintf  (fo,"%d\t%d\n",nf, 

} 

} 

} 

else  if  ((subframe  ==  -1)  «&  ( interval  !=  0))/*  Handle  super 

{  for  0=0;  j<=nwpf[fi];  j+=interval) 

{  data  =  (unsigned)out.data[word-l+j]/pltln; 
data  =  (unsigned)out.data[word-l+j]- 

plot[data]  = 
plot[data+l]  =  ’\0'; 
printf  ("%d%s\n",nf,  plot); 
plot[data]  = '  ’; 
plot[data+l]  = ' '; 

if  (result  >  2  )fprintf(fo,"%d\t%d\n",nf, 

\ 

> 

} 

else  /*  Just  plain  old  data  */ 

{  data  =  (unsigned)out.data[word-l]/pltln; 

data  =  (unsigned)out.data[word-l]-data*pltln; 
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plot[data]  = 
plot[data+l]  =  '\0'; 
printf  ("%d%s\n",nf,  plot); 
plot[data]  = ' 
plot[data+l]  = ' 

if  (result  >  2)fprintf  (fo,"%d\t%d\n",nf,  out.data[word-l]); 

} 

} 

else 

{  /*  Do  nothing  */ 

} 

printf  ("Successful  EndOfFile  reached  %d  Minor  FramesAn",  nf);/*  Cleanup,  close-up, 
and  start  over  */ 
fclose(f); 

if  (result  >  2)fclose(fo); 
nf  =  0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]); } 
printf  ("  Enter  the  file  Number  :  \n"); 
printf  ("  <ctrl>  c  to  end;  \n"); 

} 

}  /*  End  Main  */ 
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TEN  GREA 


#mclude  <stdio.h> 

#include  <string.h> 

#define  numfiles  10 

#defme  lengfiln  40 

#define  pltln  125 

int  result,  status,  subframe,  fi,  a; 

unsigned  int  data, word,  interval,  rs,  sfid; 

char  filename  [numfiles]  [lengfiln],  ofilename[IengfiIn],  ifilename[lengfiln],  outdat  [lengfiln], 
plot[pltln+l]; 

char  fsll  [numfiles]  [lengfiln]; 

unsigned  int  b=0,i,x=0,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido [numfiles], 
sff[numfiles],  sfl [numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  nfrms [numfiles],  wow[lengfiln]; 

unsigned  long  int  z=0,  j,  1,  n,  nfi=0,fsp[numfiles],  nbps[numfiles],afbl[4096],counterl[4096]; 
char  response,  excel[32]; 

char  hello,newfilel  [30],newfile2[30],newfile3[30],newfile[40]; 
struct  tbits 


unsigned  int  bl 
unsigned  int  b2 
unsigned  int  b3 
unsigned  int  b4 


struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  :  8; 


struct  kbits 
{ 


unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 


int  bl  :  2; 
int  b2  : 2; 
int  b3  :  2; 
int  b4  :  2; 
int  b5  : 2; 
int  b6  ;  2; 
int  b7  :  2; 
int  b8  :  2; 
int  b9  :  2; 
int  bl 0:2; 
int  bll:  2; 
int  b  12:  2; 
intbl3:2; 
int  bl4:  2; 
int  bl5:  2; 


unsigned  int  bl6:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsisned  int  n4  :  8; 

union 

^  struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in  ; 
union 

struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

unsigned  long  int  temp; 

FILE  *f,  *fin,  *fo,  *lpl,  *lp2,  *lp3,  *nice; 
unsigned  long  int  n=0; 
void  main() 

{  unsigned  long  int  1=0; 
i  =  0; 

fm=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpfIi],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sffIi],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin."%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n  , 

&filename  [i]  ,&fsp  [i]  ,&nwpfi;i]  ,&nbpw[i]  ,&nbps  [i]  ,&sfidw[i],&sfido[i]  ,&sffli]  ,&sfl  [i]  ,& 
fsl[i],&fsll[i],&nfnns[i])) 
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{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fm); 

for  (i=0;  i  <  numfile;  i++){  printf("%d-— >%s\n",i,filename[i]);}/*  Display  Menu  *! 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  ;  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpf[fi],nbpw[fi]); 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  -  %ld\n", 
nbps[fi],(long)  nfnns[fi]); 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sfflfi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 

newfile[0]='e'; 

newfile[l 

b=0; 

do 

^  newfile[2+b]=filename[fi][b];/*This  puts  'e;'  in  front  of  the  file 

name  so  it  can  read  it  off  a  cd*/ 

b++; 

}  while(b<=l  0);b=0; 
printf('’\nWORKING! ! !  !\n"); 

f^fopen(newfile,"rb");/*  Open  selected  data  file  */ 
n=0; 
z=0; 
do 

{ 

afbl[n]=n; 
counter!  [n]=0; 
n++; 

}  while(n<=l  024);n=0; 

while  (fread(in.data,(nwpf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1  ;/* Frame  counter*/ 
i=0; 
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for  (j=0;j<nwpftfi]*nbpw[fi]/8;j=j+5)/*  Unravel  packed  10  bit  data 


out.tbit[i].b5  =  in.tbit|j].b4; 
out.tbit[i].b4  =  in.tbit[j].b3; 
out.tbit[i].b3  =  in.tbit[j].b2; 
out.tbit[i].b2  =  in.tbit[j].bl; 
out.tbit[i].bl  =  in.tbitp+l].b4; 
out.tbit[i+l].b5  =  in.tbit[j+l].b3; 
out.tbit[i+l].b4  =  in.tbitp+l].b2; 
out.tbit[i+l].b3  =  in.tbit|j+l].bl; 
out.tbit[i+l].b2  =  in.tbit[j+2].b4; 
out.tbit[i+l].bl  =  in.tbit[j+2].b3; 
out.tbit[i+2].b5  =  in.tbit0+2].b2; 
out.tbit[i+2].b4  =  in.tbit|j+2].bl; 
out.tbit[i+2].b3  =  in.tbit[j+3].b4; 
out.tbit[i+2].b2  =  in.tbit0+3].b3; 
out.tbit[i+2].bl  =  in.tbit0+3].b2; 
out.tbit[i+3].b5  =  in.tbit[j+3].bl; 
out.tbit[i+3].b4  =  in,tbit0+4].b4; 
out.tbit[i+3].b3  =  in.tbit[j+4].b3; 
out.tbit[i+3].b2  =  in.tbit0+4].b2; 
out.tbit[i+3].bl  =  in.tbit|j+4].bl; 
i=i+4; 

} 

for(n=0;n<nwpf[fi]*nbpw[fi]/l 0;  n++) 

{ 

for(l=0;l<=1024;l=l+l) 

{ 

if(afbl  [l]=out.data[n]) 

{ 

counterl[l]++; 

break; 


printf("\n\nDone\n\n"); 

printf("What  is  the  name  of  the  output  file: "); 

scanf("%s",newfile2); 

nice=fopen(ne'wfile2,"w"); 

fprintf(nice,"%s\n",filename[fi]); 

z=0; 

do 


fprintf(nice,"%lx\tCount\t%ld\n",(long)afb  1  [z],(long)counter  1  [z]); 
z=z+l; 

}while(z<=1024);z=0; 
printf("\nFile  Saved  at  %s\n\n",newfile2); 


fclose(nice); 

fclose(f); 

fclose(lpl); 

•  fclose(lp2); 
nf  =  0; 
i=0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]);  }/*Prints  menu  of  files*/ 
printf  ("  Enter  the  file  Number  :  \n"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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TWOONE 


#include  <stdio.h> 

#include  <string.h> 

#define  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles][lengfiln],  ofilename[lengfiln],  ifilenaine[lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 

char  fsll[numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data, word,  interval,  rs,  sfid; 

unsigned  int  i,  j,  z=0,x=0,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido[numfiles], 
sff[numfiles],  sfl[numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  n=0,  nfrms [numfiles],  wow[lengfiln]; 
unsigned  long  int  nf=0,fsp [numfiles], 

nbps[numfiles],afbl[256],counterl[256],afb2[256],counter2[256],afb3[256],counter3[256]; 
char  response,  excel[32]; 

char  hello,newfiIel  [30],newfile2[30],newfile3[30]; 
struct  tbits 
{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  ;  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 


struct  bytes] 

unsigned 

}; 

struct  kbits 

{ 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 

unsigned 


int  nl  ;  8; 


int  bl  :  2; 
int  b2  :  2; 
int  b3  :  2; 
int  b4  : 2; 
int  b5  :  2; 
int  b6  ;  2; 
int  b7  :  2; 
int  b8  : 2; 
int  b9  :  2; 
int  blO:  2; 
int  bl  1:  2; 
intbl2:2; 
int  bl3:  2; 
intbl4:2; 


unsigned  int  bl 5:  2; 
unsigned  int  b  16:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in ; 

union 

struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  lb)^es  byte  [4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fin,  *fo,  *lpl,  *Ip2,  *lp3; 

void  mainQ 
{i  =  0; 

fin=fopen("cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl  [i]  ,&fsll  [i]  ,&nfrms  [i] ) ; 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n  , 

&filename  [i]  ,&fsp  [i]  ,&nwpfli]  ,&nbpw[i]  ,&nbps  [i]  ,&sfidw[i]  ,&sfido  [i]  ,&sffLi]  ,&sfl  [i]  ,& 
fsl[i],&fsll[i],&nfrms[i])) 

{  if  (i  <  numfiles-1)  i  =  i  +1 ; 
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} 

numfile  =  i ; 
fclose  (fm); 

for  (i=0;  i  <  numfile;  i++){  printf("%d— ->%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d’’,&fi);  fi  !=  -1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

■  {  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpflfi],nbpw[fi]); 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi]  ,sfido  [fi]  ,sff[fi]  ,sfl  [fi] ) ; 

printf^Number  of  Bits  Per  Second  =  %Id\tTotaI  Number  of  Frames  =  %ld\n", 
nbps[fi],(long)  nfrms[fi]); 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sfflfi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l ; 


zero*/ 


zero*/ 


printfC'HIT  RETURN  TO  GO  ON! ! !"); 
hello=getchar(); 
printf("\nWORKING!!  !\n"); 

f=fopen(&filename[fi][0],"rb");/*  Open  selected  data  file  */ 
n=0; 
z=0; 
do 

^  afb2[z]=z;/*This  puts  numbers  into  the  array  for  comparing*/ 
counter2[z]=0;/*This  puts  zeros  in  the  array  to  set  the  count  to 

Z++; 

}while(z<=15);z=0;/*Puts  numbers  up  to  15  in  the  array*/ 
do 

^  afbl  [n]=n;/*This  puts  numbers  into  the  array  for  comparing*/ 
counterl  [n]==0;/*This  puts  zeros  in  the  array  to  set  the  count  to 


n++; 

}while(n<=255);n=0;/*Puts  numbers  up  to  255  in  the  array*/ 
while  (fread(in.dataXnwpfIfi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  I  ;/*Frame  counter*/ 
j=0; 

for  (i=0;  j<nwpfIfi]*nbpw[fi]/8;  j++) 

{ 
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match*/ 


do 

{ 

if(in.dataO]==afbl  [n])/*Check  if  the  numbers 

{ 

counter!  [n]+=l  ;/*counter  if  they 


match*/ 


bit  data  into  nibles*/ 


frame  of  12  bit  words*/ 


break;/*Breaks  once  counted*/ 

} 

n-H-; 

}while(n<=255);n=0; 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  12 

{  out.nible[i]  .n3  =  in.nible[j]  .n2; 
out.nible[i]  .n2  =  in.nible|j]  .nl; 
out.nible[i]  .nl  =  in.nible[j+l].n2; 
out.nible[i+l].n3  =  in.nible[j+l].nl; 
out.nible[i+l  ]  .n2  =  in.nible[j+2]  .n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 
i=i+2; 

} 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/12;  j++)/*Loop  through  the 
{  z=0; 

do 


{ 

if(out.nible[j].nl==afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nible[j].n2=afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nible[j].n3==afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it.  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 


Z-H-; 

}while(z<=l  5);z=0; 


} 

printf("\n%s  DONE!  1  !!\n",filename[fi]); 
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printf("\nOutput  file  name  for  the  two  hex  #:  "); 
scanf("%s",newfilel); 

printf("\nOutput  file  name  for  the  one  hex  #: "); 
scanf("%s",newfile2); 

lp2=fopen(newfile2,"w");/* Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lp2,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
lpl=fopen(newfilel,''w");/*Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lpl,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
n=0; 
do 
{ 

fprintf(lpl,"%lx\tCount\t%ld\n",afbl  [n],(long)counterl  [n]);/* Printing  to  the 

output  file*/ 

n++; 

}while(n<=255);n=0; 

z=0; 

do 

{ 

fprintf(lp2,"%lx\tCount\t%ld\n",afb2[z],(long)counter2[z]);/*Printing  to  the 

output  file*/ 

Z++; 

}while(z<=15);z=0; 

printf  ("\nSuccessfuI  EndOfFile  reached  %ld  Minor  Frames.\n\n",(long)  nf);/*  Cleanup, 
close-up,  and  start  over  */ 

fclose(f); 
fclose(lpl); 
fclose(lp2); 
nf  =  0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]);  }/*Prints  menu  of  files*/ 
printf  ("  Enter  the  file  Number  :  \n"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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TWOONECD 


#include  <stdio.h> 

#include  <string.h> 

#define  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles] [lengfiln],  ofilename[lengfiln],  ifilename[lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 

char  fsll[numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data,word,  interval,  rs,  sfid; 

unsigned  int  b=0,i,  j,  z=0,x=0,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles], 
sfido[numfiles], 

sff[numfiles],  sfl [numfiles],  numfile,  fsl[numfiles],  nb; 
unsigned  int  n=0,  nfrms [numfiles],  wow[lengfiln]; 
unsigned  long  int  nf=0,fsp[numfiles], 

nbps[numfiles],afbl  [256], counter  l[256],afb2[256],counter2[256],afb3[256],counterj  [256]; 
char  response,  excel[32]; 

char  hello,newfilel[30],newfile2[30],newfile3[30],  newfile[40];, 
struct  tbits 
{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  :  8; 

}; 

struct  kbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  ;  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  blO;  2; 
unsigned  intbll:  2; 


unsigned  int  b  12:  2; 
unsigned  int  bl3:  2; 


unsigned  int  bl4:  2; 
unsigned  int  bl5:  2; 
unsigned  int  b  16:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  Ib54;es{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in ; 
union 


struct  kbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 
unsigned  long  int  data[4096]; 

}  out; 

FILE  *f,  *fin,  *fo,  *lpl,  *lp2,  *lp3; 


void  main() 

{i  =  0; 

fin=fopen("e:cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpfIi],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sffIi],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fm,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename  [i]  ,&fsp  [i]  ,&nwpf[i]  ,&nbpw[i]  ,&nbps  [i]  ,&sfidw[i]  ,&sfido  [i]  ,&sff[i]  ,&sfl  [i]  ,(&: 
fsl  [i]  ,&fsll  [i]  ,&nfrms[i])) 

101 


{  if  (i  <  numfiles-1)  i  -  i  +1 ; 

} 

numfile  =  i ;  . 
fclose  (fin); 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n’',iJilename[i]);}/*  Display  Menu  *! 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi],fsp[fi],nwpflfi],nbpw[fi]); 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],(long)  nfrms[fi]); 

printfC'Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=I; 

while  ((sfl[fi]-sff[fi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 

newfile[0]='e'; 

newfile[l]=':'; 

b=0; 

do 

^  newfile[2+b]=filename[fi][b];/*This  puts  'e:'  in  front  of  the  file 
name  so  it  Ccin  read  it  off  a  cd*/ 

b++; 

}while(b<=10);b=0; 
printf("\n  WORKING!  I !  !\n"); 

f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 
n=0; 
z=0; 
do 

afb2[z]=z;/*This  puts  numbers  into  the  array  for  comparing*/ 
counter2[z]==0;/*This  puts  zeros  in  the  array  to  set  the  count  to 

zero*/ 

Z++; 

}while(z<=15);z=0;/*Puts  numbers  up  to  15  in  the  array*/ 
do 

afbl  [n]=n;/*This  puts  numbers  into  the  array  for  comparing*/ 
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zero*/ 


match*/ 


counter  1  [n]=0;/*This  puts  zeros  in  the  array  to  set  the  count  to 

n++; 

}while(n<=255);n=0;/*Puts  numbers  up  to  255  in  the  array*/ 
while  (fread(in.data,(nwpflfi]*nbpw[fi]/8),Lf)y*  Read  data  records  until  EOF  */ 
{  nf  =  nf  +  1  ;/*Frame  counter*/ 
j=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j++) 

{ 

do 

{ 

if(in.dataO]==afbl[n])/*Check  if  the  numbers 

{ 

counter  l[n]+=l;/*  counter  if  they 


match*/ 


bit  data  into  nibles*/ 


break;/*Breaks  once  counted*/ 

} 

n++; 

}while(n<=255);n=0; 

i=0; 

for  0=0;  j<nwpfi;fi]*nbpw[fi]/8;j+=3)/*  Unravel  packed  12 

{  out.nible[i]  .n3  =  in.nible|j]  .n2; 
out.nible[i]  .n2  =  in.nible[j]  .nl; 
out.nible[i]  .nl  =  in.nible0+l].n2; 
out.nible[i+l].n3  =  in.nible[j+l].nl; 
out.nible[i+l].n2  =  in.nible[]+2].n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 
i=i+2; 


frame  of  12  bit  words*/ 


} 

for  0=0;  j<nwpf[fi]*nbpw[fi]/12;  j++)/*Loop  through  the 
{  z=0; 

do 


{ 

if(out.nibleO].nl==afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nible[i].n2==afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 

{ 

counter2[z]+=l; 

} 

if(out.nible[j].n3==afb2[z])/*Checks  if  the 
numbers  match  up  and  then  counts  it,  goes  bit  by  bit*/ 
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{ 

counter2[z]+=l; 

} 

Z++; 

}  while(z<=  15)  ;z=0 ; 


} 

printf("\n%s  DONE! ! I !\n",filename[fi]); 
printf("\nOutput  file  name  for  the  two  hex  #: "); 
scanf("%s",newfile  1 ); 

printf("\nOutput  file  name  for  the  one  hex  #: "); 
scanf("%s",newfile2); 

lp2=fopen(newfile2,"w");/*Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lp2,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
lpl=fopen(newfilel,"w");/*Opens  a  file  to  put  the  numbers  and  the  count  into*/ 
fprintf(lpl,"%s\n",filename[fi]);/*Prints  file  name  to  the  output  file*/ 
n=0; 
do 

^  fprintf(lpl,’'%lx\tCount\t%ld\n",afbl[n],(long)counterl[n]);/*Printing  to  the 

output  file*/ 

n++; 

}while(n<=255);n=0; 

z=0; 

do 

^  fprintf(lp2,"%Ix\tCount\t%ld\n",afb2[z],(long)counter2[z]);/*Printing  to  the 

output  file*/ 

Z++; 

}while(z<=l  5);z=0; 

printf  ("\nSuccessful  EndOfFile  reached  %ld  Minor  Frames.\n\n",(long)  nf);/*  Cleanup, 
close-up,  and  start  over  */ 

fclose(f); 

fclose(lpl); 

fclose0p2); 

nf=0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]);  }/*Prints  menu  of  files*/ 
printf  ("  Enter  the  file  Number  :  \n"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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WORD 


/* 


* 

* 

Name: 

Display_Data.c 

* 

Creation  Date: 

09-Sep-97 

* 

Description: 

Display  Data  from  Telemetry  Sample  Data  Sets 

* 

Author: 

Vernon  Diekmann,  TYBRIN 

* 

Target  Machine/OS/Compiler: 

Intel  Pentium/Windows  95/InT/MS  Visual  C++ 

* 

* 

Revision  History 

♦ 

Initials  Date 

Description 

* 

VD 

09/09/97 

Initial  version 

* 

VD 

10/10/97 

Added  12  and  16  bit 

* 

VD 

01/15/98 

Added  24  bit  and  read  byte  records 

* 

* 

VD 

01/28/98 

Input  the  file  name  for  the  data.dat  file. 

*/ 

#include  <stdio.h> 

#include  <string.h> 

#define  numfiles  10 
#define  lengfiln  40 
#define  pltln  125 

char  fllename[numfiles] [lengfiln],  ofilename [lengfiln],  ifilename[lengfiln],  outdat[lengfiln], 
plot[pltln+l]; 

char  fsll  [numfiles]  [lengfiln]; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data, word,  interval,  rs,  sfid; 

unsigned  int  i,  j,  x=0,nwpflnumfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido [numfiles], 
sff[numfiles],  sfl[numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  n=0,  nfrms[numfiles],  wow[lengfiln]; 
unsigned  long  int  z=0,nf=0,fsp[numfiles], 

nbpslnumfiles],afbl[256],counterl[256],afb2[4096],counter2[4096],afb3[256],counter3[256]; 
char  response,  excel[32]; 

char  hello,newfile  1  [3  0]  ,newfile2  [3  0]  ,newfile3  [30]; 
struct  tbits 
{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
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unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  ;  8; 

struct  kbits 

{ 

unsigned  int  bl  ;  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  ;  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  :  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  b  10:  2; 
unsigned  int  bl  1:  2; 
unsigned  int  bl2:  2;  ’ 
unsigned  int  bl3:  2; 
unsigned  int  b  14:  2; 
unsigned  int  bl5:  2; 
unsigned  int  bl6:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  lbytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  16384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in ; 

union 
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{ 


struct  Itbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  -  4096*/ 

struct  Inibles  nible[4096]; 
struct  Ibytes  byte[4096]; 
unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fin,  *fo,  *lpl,  *lp2,  *lp3; 
void  mainO 


for  (i  =  0;  i  <  pltln;  i++)  pIot[i]  = '  ';/*  Initialize  array  */ 
i  =  0; 

fin=fopen("cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sffIi],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i])) 

{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fin); 

for  (i==0;  i  <  numfile;  i++){  printf("%d- — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  :  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",&fi);  fi  !=  -1;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi]  ,fsp  [fi]  ,nwpf[fi]  ,nbp  w[fi] ) ; 

printf("Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi]  ,sfido[fi]  ,sff[fi]  ,sfl  [fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],(long)  nfrms[fi]); 

printf("Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sfflfi])>i)  {nb=nb+l;i=2*i;} 
rs  =  nbpw[fi]-sfido[fi]-nb+l; 
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printf  ("  Enter\nF  to  print  each  frame  or  \nW  to  plot  a  single  word:  \n"); 
response  =  getchar(); 


if  (response  =  'w'  1|  response  =  'W') 

printf  ("  Enter  the\n  Word,  Subframe,  Interval  to  display:  \n"); 
result  =  scanf("%d%d%d",&word,&subframe,&interval); 
hello  =  getcharO; 

printf  ("  Enter  'c'  to  continue  or  type  the  path  and  file  name  for  the  data  of 
the  word  selected:  \n"); 

result  =  scanf("%s",&excel); 

} 

f=fopen(&filename[fi][0],"rb");/*  Open  selected  data  file  */ 
result  =  strlen(&excel[0]); 
if(result  >  2) 

fo=fopen(&excel[0],"w");/*  Open  option  data  output  file  for  excel  */ 
j=0; 
z=0; 
do 
{ 

afb2[z]=z; 

counter2[z]=0; 

Z++; 

}while(z<=4095);z=0; 

printf("\nHello\n"); 

while  (fread(in.data,(nwpf[fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 

{  nf  =  nf  +  1;  printf("%ld^",nf); 

if  (response  ==  'f  ||  response  'F')/*  Display  each  minor  frame  */ 

(  /*  *printf  ("\nFrame=%d\n",nf);  *  */ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  12 

bit  data  */ 

{  out.nible[i]  .n3  =  in.nible[j]  .n2; 
out.nible[i]  .n2  =  in.nible[j]  .nl; 
out.nible[i]  .nl  =  in.nible0+l].n2; 
out.nible[i+l].n3  =  in.nible[j+l].nl; 
out.nible[i+l].n2  =  in.nible[j+2].n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 
i=i+2; 

} 

for(i=0;  i<nv^f[fi]*nbpw[fi]/12;  i++) 

{  z=0; 

do 

{ 


if(out.data[i]==afb2[z]) 
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parameter  */ 


bit  data  */ 


bit  data  */ 


{ 

counter2[z]+=l; 

break; 

} 

Z-H-; 

}while(z<=4095);z=0; 

} 

} 

else  if  (response  ==  'w'  |1  response  ==  'W')  /*  Display  selected  data 

{ 

if  (nbpw[fi]==8) 

{ 

i=0; 

for  (j=0;  j<nwpftfi]*nbpw[fi]/8;  j+=l)/*  Unravel  packed  8 


{ 

out.byte[i].nl  =  in.byte[j].nl; 

} 

} 

if  (nbpw[fi]==10) 

{ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=5)/*  Unravel  packed  .10 


{ 

out.tbit[i].b5  =  in.tbit(j]-b4; 

out.tbit[i].b4  =  in.tbit[j].b3; 

out,tbit[i].b3  =  in.tbit[}].b2; 

oul.tbit[i].b2  =  in.tbit[j].bl; 

cui.tbit[i].bl  =  in.tbit0+l].b4; 

out.tbit[i+l].b5  =  in.tbit[j+l].b3; 
out.tbit[i+l].b4  =  in.tbit[j+l].b2; 
out.tbit[i+l].b3  =  in.tbit|j+l].bl; 
out.tbit[i+l].b2  =  in.tbit0+2].b4; 
out.tbit[i+l].bl  =  in.tbit|j+2].b3; 
out.tbit[i+2].b5  =  in.tbit0+2].b2; 
out.tbit[i+2].b4  =  in.tbit[j+2].bl; 
out.tbit[i+2].b3  =  in.tbit[j+3].b4; 
out.tbit[i+2].b2  =  in.tbit[j+3].b3; 
out.tbit[i+2].bl  =  in.tbit[j+3].b2; 
out.tbit[i+3].b5  =  in.tbit[)+3].bl; 
out.tbit[i+3].b4  =  in.tbit[j+4].b4; 
out.tbit[i+3].b3  =  in.tbit[j+4].b3; 
out.tbit[i+3].b2  =  in.tbit[j+4].b2; 
out.tbit[i+3].bl  =  in.tbit0+4].bl; 
i=i+4; 
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} 


if  (nbpw[fi]==12) 

{ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  12 

bit  data  */ 

{  out.nible[i]  .n3  =  in.nibleU]  .n2; 
out.nible[i]  .n2  =  in.nible[j]  .nl; 
out.nible[i]  .nl  =  in.nible[jM].n2; 
out.nible[i+l].n3  =  in.niblefj  HJ.nl; 
out.nible[i+l].n2  =  in.nible[i+2].n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 
i=i+2; 

} 

} 

if  (nbpw[fi]==16) 

{ 

i=0; 

for  (j=0;  j<nwpf{fi]*nbpw[fi]/8;  j+=2)/*  Unravel  packed  16 

bit  data  */ 

{  out.byte[i].n2  =  in.byte[j]  .nl; 

out.byte[i].nl  =in.byte[j+l].nl; 
i=i+l; 

} 

} 

if  (nbpw[fi]==24) 

{ 

i=0; 

for  (j=0;  j<nwpf[fi]*nbpw[fi]/8;  j+=3)/*  Unravel  packed  24 

bit  data  */ 

{  out.byte[i].n3  =  in.byte|j]  .nl; 

out.byte[i].n2  =  in.byte|j+l].nl; 
out.byte[i].nl  =  in.byte[j+2].nl; 
i=i+l; 

} 

} 

if  ((sfidw[fi]  >=  0)  &  ( interval  !=  0)  &  (subframe  >=  0))/*  Handle 

sub  frame  data  */ 

{  sfid  =  out.data[sfidw[fi]-l  ]; 

sfid  =  sfid  »  rs; 

for  (j=0;  j<(sfl[fi]-sff[fi]+l);  j+=interval) 

{  if  ((subframe+j)  ==  sfid) 

{  data  =  (unsigned)out.data[word-l]/pltln; 
data  =  (unsigned)out.data[word-l]- 

data*pltln; 

plot[data]  = 
plot[data+l]  =  '\0'; 
printf  ("%d%s\n",nf,  plot); 

110 


plot[data]  = ' 
plot[data+l]  = ' 

if  (result  >  2)fprintf  (fo,"%d\t%d\n",nf, 


out.data[word-l]); 


commutated  data  */ 


data*pltln; 


out.data[word-l+j]); 


} 

} 

} 

else  if  ((subframe  =  -!)&(  interval !-  0))/*  Handle  super 

{  for  (j=0;  j<=nwpf[fi] ;  j+=interval) 

{  data  =  (unsigned)out.data[word-l+j]/pltln; 
data  =  (unsigned)out.data[word-l+j]- 

plot[data]  = 
plot[data+l]  =  ’\0'; 
printf  ("%d%s\n",nf,  plot); 
plot[data]  = ' 
plot[data+l]  = ' 

if  (result  >  2  )fprintf(fo,"%d\t%d\n",nf. 


} 


} 


} 

else 

{ 

} 


} 

else  /*  Just  plain  old  data  */ 

{  data  =  (unsigned)out.data[word- 1  J/pltln; 

data  =  (imsigned)out.data[word-l]-data*pltln; 

plot  [data]  = 

plot[data+l]  =  '\0'; 

printf  ("%d%s\n",nf,  plot); 

plot[data]  =  "; 

plot[data+l]  =  "; 

if  (result  >  2)fprintf  (fo,"%d\t%d\n",nf,  out.data[word-l]); 


/*  Do  nothing  */ 


printf("\nOutput  file  name  for  the  word  data; "); 
scanf("%s’',newfile2); 
lp2=fopen(newfile2 , "  w"); 
fprintf(lp2,"%s\n",filename[fi]); 

z=0; 

do 

{ 
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fyrintf(lp2,"%lx\tCount\t%ld\n",afb2[z],(long)  counter2[z]); 

Z++; 

}while(z<=4095);z=0; 

printf  ("Successful  EndOfFile  reached  %ld  Minor  Frames.\n",  (long)nf);/*  Cleanup, 
close-up,  and  start  over  */ 
fclose(f); 
fclose(lp2); 

if  (result  >  2)fclose(fo); 
nf=0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]); } 
printf  ("  Enter  the  file  Number  :  \n"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main*/ 
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#include  <stdio.h> 

#include  <string.h> 

#defme  numfiles  1 0 
#define  lengfiln  40 
#define  pltln  125 

char  filename[numfiles][lengfiln],  ofilename[lengfiln],  ifilename[lengfiln],  outdat[lengfiln], 
plot[pltIn+l]; 

char  fsll  [numfiles]  [lengfilnj; 
int  result,  status,  subframe,  fi,  a; 
unsigned  int  data,word,  interval,  rs,  sfid; 

unsigned  int  i,  j,  x=0,nwpf[numfiles],  nbpw[numfiles],  sfidw[numfiles],  sfido [numfiles], 
sff[numfiles],  sfl [numfiles],  numfile,  fsl [numfiles],  nb; 
unsigned  int  b=0,n=0,  nfrms [numfiles],  wow[lengfiIn]; 
unsigned  long  int  z=0,nf=0,fsp[numfiles], 

nbpsl;numfiles],afbl[256],counterl[256],afb2[4096],counter2[4096],afb3[256],counter3[256] 

char  response,  excel[32]; 

char  hello,newfilel[30],newfile2[30],ne\vfile3[30],  newfile[40]; 
struct  tbits 


{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  :  2; 
unsigned  int  b4  :  2; 

}; 

struct  nibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 

}; 

struct  bytes  { 

unsigned  int  nl  :  8; 

}; 

struct  kbits 

{ 

unsigned  int  bl  :  2; 
unsigned  int  b2  :  2; 
unsigned  int  b3  ;  2; 
unsigned  int  b4  :  2; 
unsigned  int  b5  :  2; 
unsigned  int  b6  :  2; 
unsigned  int  b7  ;  2; 
unsigned  int  b8  :  2; 
unsigned  int  b9  :  2; 
unsigned  int  blO:  2; 
unsigned  int  bll:  2; 
unsigned  int  b  12:  2; 
unsigned  int  bl3:  2; 
unsigned  int  b  14:  2; 


unsigned  int  bl5:  2; 
unsigned  int  bl6:  2; 

}; 

struct  lnibles{ 

unsigned  int  nl  :  4; 
unsigned  int  n2  :  4; 
unsigned  int  n3  :  4; 
unsigned  int  n4  :  4; 
unsigned  int  n5  :  4; 
unsigned  int  n6  :  4; 
unsigned  int  n7  :  4; 
unsigned  int  n8  :  4; 

}; 

struct  Ibytes{ 

unsigned  int  nl  :  8; 
unsigned  int  n2  :  8; 
unsigned  int  n3  :  8; 
unsigned  int  n4  :  8; 

}; 

union 

struct  tbits  tbit[2048];/*  Max  allowable  Class  II  bytes  per  frame  1 6384/8  = 

2048*/ 

struct  nibles  nible[2048]; 
struct  bytes  byte[2048]; 
unsigned  char  data[2048]; 

}  in ; 
union 


struct  Itbits  tbit[4096];/*  Max  allowable  Class  II  wpf  =  16384/4  =  4096*/ 
struct  Inibles  nible[4096]; 
struct  Ibytes  byte  [4096]; 

unsigned  long  int  data[4096]; 

}  out ; 

FILE  *f,  *fin,  *fo,  *lpl,  *lp2,  *lp3; 


void  mainQ 
{  i  =  0; 

fin=fopen("e;cinput.dat","rb");  /*  Open  file  and  read  data  description  information  */ 
for  (status  =  fscanf(fln,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 

&filename[i],&fsp[i],&nwpf[i],&nbpw[i],&nbps[i],&sfidw[i],&sfido[i],&sff[i],&sfl[i],& 

fsl[i],&fsll[i],&nfrms[i]); 

status  !=  EOF; 

status  =  fscanf(fin,"%s  %ld  %d  %d  %ld  %d  %d  %d  %d  %d  %s  %ld\n", 


&filename  [i]  ,&fsp  [i]  ,&nwpf[i]  ,&nbpw[i]  ,&nbps  [i]  ,&sfidw[i],&sfido[i]  ,&sffli]  ,&sfl  [i]  ,& 
fsl[i],&fsll[i],&nfnns[i])) 
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{  if  (i  <  numfiles-1)  i  =  i +1; 

} 

numfile  =  i ; 
fclose  (fin); 

for  (i=0;  i  <  numfile;  i++){  printf("%d — >%s\n",i,filename[i]);}/*  Display  Menu  */ 
i=0;  fi  =  0;  result  =  0; 

printf  ("  Enter  the  file  Number  ;  \n");/*  Make  Selection  */ 

for  (result=scanf("%d",«S:fi);  fi  !=  -1 ;  result=scanf("%d",&fi))/*  Main  loop  executed 
once  for  each  file  selection  */ 

{  hello  =  getcharO;  nf  =  0;/*  First  display  and  collect  option  information  */ 

printf("\nFile  =  %s\t\tFrame  Sync  Pattern  =  %lx\nNumber  of  Words  Per  Frame  = 
%d\t\tNumber  of  Bits  Per  Word  =  %d\n", 

filename[fi]  ,fsp[fi]  ,nwpf[fi]  ,nbpw[fi] ) ; 

printfC'Sub  Frame  ID  Word  =  %d\t\t\tSub  Frame  ID  Offset  =  %d\nSub  Frame 
First  =  %d\t\t\tSub  Frame  Last  =  %d\n", 

sfidw[fi],sfido[fi],sff[fi],sfl[fi]); 

printf("Number  of  Bits  Per  Second  =  %ld\tTotal  Number  of  Frames  =  %ld\n", 
nbps[fi],(long)  nfrms[fi]); 

printfC’Frame  Sync  Length  =  %d\t\tSync  Location  =  %s\n\n", 

fsl[fi],fsll[fi]); 

nb=0;i=l; 

while  ((sfl[fi]-sfflfi])>i)  {nb=nb+l;i=2*i;} 

rs  =  nbpw[fi]-sfido[fi]-nb+l; 

newfile[0]='e'; 

newfile[l]=':'; 

b=0; 

do 

^  newfile[2+b]=filename[fi][b];/*This  puts  'e;'  in  front  of  the  file  so 

it  can  read  off  a  cd*/ 

b++; 

}  while(b<=  1 0);b=0 ; 

printf("\nWORKING! ! !  !\n"); 
f=fopen(newfile,"rb");/*  Open  selected  data  file  */ 

j=0; 

z=0; 

do 

{ 

afb2[z]=z; 

counter2[z]=0; 

Z++; 

}while(z<=4095);z=0; 

while  (fread(in.data,(nwpfi;fi]*nbpw[fi]/8),l,f))/*  Read  data  records  until  EOF  */ 

{  nf=nf+l; 

i=0; 

for  G=0;j<nwpf[fi]*nbpw[fi]/8;j+=3)/*  Unravel  packed  12 

bit  data  */ 


115 


{  out.nible[i]  .n3  =  in.nible[j]  .n2; 
out.nible[i]  .n2  =  in.nible[j]  .nl; 
out.nible[i]  .nl  =  in.nible[j+l].n2; 
out.nible[i+l].n3  =  in.nible[j+l].nl; 
out.nible[i+l].n2  =  in.nible[i+2].n2; 
out.nible[i+l].nl  =  in.nible[j+2].nl; 
i=i+2; 

} 

for(i=0;  i<nwpf[fi]*nbpw[fi]/12;  i++) 

{  z=0; 

do 

{ 

if(out.data[i]==afb2  [z]) 

{ 

counter2[z]+=l; 


break; 

} 


Z++; 

}while(z<=4095);z=0; 

} 

} 

printf("\n\nFILE  %s  DONE\n\n",filename[fi]); 

printf("\nOutput  path  and  file  name  for  the  word  data: "); 

scanf("%s",newfile2); 

lp2=fopen(newfile2,"w"); 

fprintf(lp2 , "  %s\n " ,  filename  [fi] ) ; 


z=0; 

do 

{ 

fprintf(lp2,"%lx\tCount\t%ld\n",afb2[z],(long)  counter2[z]); 

Z++; 

}while(z<=4095);z=0; 

printf  ("Successful  EndOfFile  reached  %ld  Minor  Frames.\n",  (long)nf);/*  Cleanup, 
close-up,  and  start  over  */ 
fclose(f); 
fclose(lp2); 

if  (result  >  2)fclose(fo); 
nf  =  0; 

for  (i=0;  i  <  numfile;  i++){  printf("%d  %s\n",i,filename[i]); } 
printf  ("  Enter  the  file  Number  :  hi"); 
printf  ("  <ctrl>  c  to  end:  \n"); 

} 

}  /*  End  Main  */ 
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Shakti  Davis  wrote  code  to  implement  the  channel  simulations  using  the  error  files  supplied  by 
TYBRIN.  The  compression  algorithm  (Izss  is  from  the  Mark  Nelson  data  compression  book). 
These  computer  programs  and  instructions  are  listed  below: 

Documentation  for  BERA  Simulations: 

To  run  a  BERA  simulation,  open  MATLAB  and  type  “bera_sim”  at  the  command 
prompt.  This  should  bring  up  a  dialog  box  in  which  you  can  specify  the  transmit  data  file,  the 
receive  data  file  and  whether  you  want  the  data  to  be  compressed  (with  Izss)  before  the  errors  are 
injected.  -The  last  parameter,  channel  coding  isn’t  implemented  yet.  To  get  the  simulation 
slcirted,  you  can  either  hit  OK  to  accept  the  default  parameters,  or  you  can  change  any  of  the  first 
three  parameters.  For  the  filenames,  I  would  recommend  giving  the  full  path  name  (e.g. 
g:\sdsdata\sds008c.dat  instead  of  just  sds008c.dat)  to  insure  that  the  simulation  uses  the  correct 
file,  Notice  that  the  files  can  be  read  from  the  CD  ROM  drive,  so  you  won’t  have  to  copy  the 
files  to  the  harddrive.  Please  be  patient  with  the  simulation  because  it  may  take  several  minutes 
to  complete.  If,  however,  the  simulation  is  running  and  no  harddrive  activity  occurs  for  several 
minutes,  then  the  program  may  have  frozen.  In  this  case,  I  would  suggest  hitting  CTRL-ALT- 
DEL  and  ending  the  MATLAB  session. 

The  simulation  is  composed  of  several  files  acting  on  the  data  files.  First,  the 
“BERA_sim.m”  MATLAB  file  collects  the  simulation  parameters  from  the  user  and  determines 
whether  the  data  file  should  be  compressed.  If  the  user  indicated  “yes”  for  compression,  then  the 
C-file  Izss.c  is  used  to  compress  the  file.  This  is  possible  by  compiling  the  C  file  as  a  MEX 
function.  To  compile  the  Izss.c  file  as  a  MEX  file,  do  the  following: 

1 .  Open  MATLAB 

2.  From  the  prompt,  change  the  working  directory  path  to  D:\stats\simulink 

simulationMzss  which  is  where  the  C  compression  code  resides 

»  cd  'D;\stats\simulink  simulationMzss' 

3.  Compile  the  compression  routine  as  a  MEX  file: 

»mex  -DMATLAB  -output' ‘d:\stats\simulink  simulationMzss_c’  main-c.c  Izss.c 
■  bitio.c  errhand.c 

this  creates  a  file  called  lzss_c.dll  in  the  d:\stats\simulink  simulation\  directory  that  can  be 

called  from  MATLAB 

4.  Use  this  function  in  a  MATLAB  script  as  a  normal  function  call 

»lzss_c('SDS001c.dat’,'SDS001c.cmp'); 

So  the  BERA_sim  file  compiles  the  input  data  file  using  this  lzss_c  MEX  function  and  moves  on 
to  the  next  step  which  is  to  process  the  BERA  file.  The  processing  is  also  done  in  a  C  program 
that  was  compiled  as  a  MEX  function.  The  C  program  is  called  procBERA.c  and  is  located  at 
d:\stats\simulink  simulation\.  To  compile  this  program  as  a  MEX  file,  do  the  following; 

1 .  Open  MATLAB 

2.  From  the  prompt,  change  the  working  directory  path  to  D:\stats\simulink  simulation 

which  is  where  the  procBERA  C  source  code  resides 

»  cd  'D:\stats\simulink  simulation\lzss' 

3.  Compile  the  compression  routine  as  a  MEX  file: 

»mex  -DMATLAB  -output  ‘d;\stats\simulink  simulation\procbera’  procbera.c 
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this  creates  a  file  called  procbera.dll  in  the  d:\stats\simulink  simulation\  directory  that  can 
be  called  from  MATLAB 

4.  Use  this  function  in  a  MATLAB  script  as  a  normal  function  call 
»procBERA('d:\stats\flight2bera\F02f_02_Ml  2_028.bin', 

'd:\stats\SDS001c.dat',  'd:\stats\rcvddata.bin') 

This  procBERA.dll  function  is  called  whether  the  compression  is  done  or  not,  so  this  MEX 
function  must  exist  in  order  to  run  the  BERA  simulation.  The  procBERA  C  code  takes  as  inputs 
the  name  of  the  BERA  error  file,  the  name  of  the  data  input  file,  and  the  name  of  the  data  output 
file.  Each  of  these  files  is  opened  and  a  replica  of  the  data  input  file  is  copied  to  the  data  output 
file  (this  is  the  easiest  way  since  there  are  so  few  errors  relative  to  the  number  of  correct  bits). 

So  originally  the  input  data  is  equal  to  the  output  data.  Next,  the  error  file  is  read  one  entry  at  a 
time,  where  each  error  entry  consists  of  the  following  6-byte  structure: 

XOR  word  :  bits  47-32 
Event  Byte  :  bits  31-24 
24  bit  counter :  bits  23-0 

For  each  error,  an  absolute  position  of  that  error  is  specified  by  the  24-bit  counter,  so  the 
program  goes  to  that  position  in  the  output  data  file.  The  XOR  word  contains  a  1  at  each  error 
location,  so  once  the  output  data  file  is  correctly  positioned,  its  data  at  that  position  (2-bytes),  can 
be  read,  XOR-ed  with  the  XOR  word  of  the  error  entry,  and  written  back  to  the  data  output  file. 
This  process  is  repeated  for  each  error  event  in  the  BERA  error  file,  or  until  the  end  of  the  data 
file  is  reached.  The  BERA  files  cover  such  a  long  transmit  signal  that  the  input  data  file  is 
usually  the  limiting  factor  (e.g.  the  position  specified  by  the  error  events  are  larger  than  the  last 
position  in  the  input  data  file).  The  event  byte  of  the  error  structure  is  not  used.  While  the  error 
file  is  being  processed,  statistics  about  the  file  are  collected.  The  three  statistics  which  are 
returned  to  the  MATLAB  program  are  the  bit-error-rate,  the  number  of  bits  in  error,  and  the  total 
number  of  bits  processed.  To  access  these  return  values  from  the  MATLAB  call,  simply  call 
procBERA  as  follows: 

[BER, numErrs, numBits] =procBERA (BERAf ile, inf ile , outf ile ) ; 
and  the  bit  error  rate  will  be  returned  in  the  BER  variable,  the  number  of  bits  in  error  will  be 
returned  in  the  numErrs  variable,  and  the  total  number  of  bits  processed  will  bw  returned  in  the 
numBits  variable.  If  the  simulation  was  specified  to  use  the  compression,  then  the  data  file  that  is 
created  by  the  procBERA  function  must  be  decompressed  using  the  lzss_e  MEX  function 
(compiled  the  same  way  as  the  lzss_c  MEX  function  but  instead  of  main-c.c  use  main-e.c).  After 
the  file  is  decompressed,  it  can  be  compared  with  the  original  data  file  and  the  true  bit-error  rate 
can  be  determined. 
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BERA  simulation  program: 


luui.  Cl  L  -L  L-ii 


usina  HERA  ril6S 


%  prompi:  for  data  file  names 

Prompts  =  {  *  Enter  the  EERA  filename:’,... 

’Enter  the  filenvame  of  the  cripinai  oata:  . 

'Enter  the  filename  for  the  output  data:',... 


o 

Q 

you  want  to  compress  the 

data  first?  (Yes/N 

'  Do 

veu  want  to  use  channel 

ceding?  (Yes/No)'}; 

{  'e 

:  \stats\fiight2bera\F02f_ 

02__M12_O23  .bin  '  ,  .  .  . 

'  e : 

\ 3 1  a  t  s \ S  DS  0  0 1 c . dat ' ,  . . . 

'  e : 

\stats\rcvddata . dat ' ,  . . . 

'No 

I 

,  .  .  . 

'No 

lineNo=l; 

title  =  'Simulation  Configuration’; 

Answer  =  inputdlg ( Prompts , title, lineNo,def); 


if  (strcmpi (char (Answer (4 )),' yes ') ) 

^  make  the  compressed  data  filename  equal  to  the  original 
V  data  filename  with  the  extension  ’'.cmp"  and  save  it  to  t: 
1  stats  directory 

index=f indstr (char (Answer (2)  ; 

if  (index) 

infile=char (Answer (2)); 
inf ile=in file (1 : index- 1)  / 
inf ile=strcat (infile, ’ .cmp' ) ; 
else 

infile=char (Answer (2) ) ; 

inf ile=strcat  (infile,  '  .  crap '  )  ; 

end 

index=findstr (char (Answer (2) ) ,  ’ \ ; 
if  (index) 

infile=strcat ( ’ e : XstatsX ’ , infile (index (end) +1 : end) ) ; 
else 

inf ile=st rcat ( ’ e : \s tats \ ' , infile ) ; 

end 


o 


l2ss_c (char (Answer (2) ) , infile) ; 

-  make  the  compressed  output  filename  equal  to  the  output 
::  filename  with  the  extension  ".cmp” 
index=f indstr (char (Answer (3) ),'.')  / 
if  (index) 

outfile=char (Answer (3)  )  ; 
outf ile=outf ile ( 1 : index- 1 )  ; 
out file=s treat (outf ile,  '  . cmp ' ) ; 
else 

outfile=char (Answer (3) ) ; 

out f ile=strcat (outf ile,  ’ . cmp ' )  / 

end 

else 

infile  =  char (Answer ( 2 )) ; 
outfile  =  char (Answer (3) ); 
end; 
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p;iroc0S5  02rjrc2r 

[BER, actualErrs, actualBits] =procBERA (char (Answer (1) ) , infile, outfile) ; 

if  (strcmpi (char (Answer ( 4 ) ) , ' yes ' ) ) 

'r decoiriDress  the  output  file  and  calculate  tne  Bii\ 
lzss_e (outfile, char (Answer (3) ) )  ; 

[BERTactualErrs, actualBits ] =checkBER ( char (Answer ( 1 ) )  ,  char (Answer (2 ) )  ) 

end 

fprintf ( 1 ,  ’ Scenario : \n ' )  ; 

fprintf  (1,  '  \tlnp''i-.  file:  \t\t  \t%s\n' ,  char  (Answer  (1)  )  )  ; 
fprintf (1, '\tOutput  f ile : \ t \L\t%3\n char (Answer (2 ) ) ) ; 
fprintf  (1,  '\tData  Coir:pressed?\t\t%3\n’ ,  char  (Answer  (4  )  )  )  ; 
fprintf  (1,  '  \tChannel  Coding?\t\t\t 'issXn ' ,  char  (Answer  (5)  )  )  ; 
fprintf (1, 'Statistics : \n' ) ; 

fprintf (1, ' \tBit  Error  Rate;\t\t\tie\n',BER); 

fprintf  (1,  '  itNumber  of  Bits  Proce;?sed: \t'sd\n' ,  actualBits)  ; 

fprintf  ( 1 ,  '  XtNumber  of  Errors  :  \t  \t  :-d\r. ' ,  actualErrs )  ; 


PROCBERA  program  (the  error  file  was  supplied  by  SyntheSys  Research); 


ifndef  _B  A25READ_H_ 

#defme  _BA25READ_H_ 

/* 

* _ * 

*  Global  Constant  Definitions 

*  _ * 

*! 

#define  ROLL_PERIOD_25  0x01000000 

#define  ROLL_BYTES_25  (ROLL_PERIOD_25  «  1) 

#define  ROLL_M ASK_25  (ROLL_B YTES_25  - 1 ) 

#define  WORD_SIZE_25  16 
I* 

^ _ * 

*  Global  Macro  Definitions 

^  _ _ _ _ _ _ _ ^ 

*1 

#define  HW25_EVENT_BYTE(Ptr)  (*(((unsigned  char  *)(Ptr))+3)) 
/* 

_ _ _ _ _ 

*  Type  Definitions 

sje  _ _ _ _ _ _ --—-I _ _  ^ 

*1 

/* 

*  BitAlyzer  25  Event  Package  Structure 

* 

*  Note:  Xor  is  declared  as  a  short  instead  of  a  bit  field 

*  so  that  sizeof(Sample25)  is  correctly  evaluated  to 

*  6  bytes. 

*/ 


typedef  struct  { 

unsigned  long  Clock  ;24; 
unsigned  long  Marker  :  1 ; 
unsigned  long  Pattern  :  1 ; 
unsigned  long  Resync  :  1 ; 
unsigned  long  Squelch  ;  1 ; 
unsigned  long  MarkerPosition  :4; 
unsigned  short  Xor ; 

}  Sample25  ; 

typedef  struct  { 

unsigned  long  data  ;16; 

}Word 


#endif 
/*  EOF  */ 


/*  BA25READ.C  */ 
/* 


===*/ 

/*  The  objective  of  this  program  is  to  show  how  a  BitAlyzer  25  Error  File 
could  be  read  using  a  simple  program  written  in  C.  To  understand  this^ 
program  the  user  must  be  familiar  with  the  BitAlyzer  error  file  forrnat. 

This  program  reads  an  error  file  and  computes  the  total  number  of  errors 
and  the  total  error  free  interval. 

For  further  information  on  writing  programs  to  access  error  file  data 
or  error  file  format  specifications  please  contact:  SyntheSys  Research, 

Inc.  at  650-364-1853 

Modified  11/18/99 -S.D. 

Now  this  program  locates  the  errors  and  actually  XORs  the  error  bits 
with  the  true  data.  The  resulting  file  can  then  be  compared  to  the 
original  data  file  for  error  analysis. 

Syntax:  procBERA  BERAfile  binary  datafile  outputfile 

Notes  from  BYU: 


The  BitAlyzer  stores  it's  data  in  6  byte  chunks. 

XOR  word  :  bits  47-32 
Event  Byte  :  bits  31-24 
24  bit  counter  :  bits  23-0 

Because  the  data  was  written  by  a  PC,  and  we  were  doing  our  work  on  Unix 
boxes,  we  had  a  little-endian  big-endian  problem  (i.e.  the  byte 
addressing  in  the  file  is  opposite  between  the  PC  and  Unix  boxes.)  To 
overcome  this  problem  I  just  read  each  byte  out  individually,  and 
processed  them  the  way  I  wanted.  I  hope  that  the  documentation  in  the 
code  is  clear.  Another  reason  to  open  the  data  in  Byte  chucks  is  since 
the  data  is  6  Bytes  long,  my  attempts  to  try  and  read  a  6  Byte  block 
were  actually  executed  as  a  read  of  6  Bytes,  then  skip  to  the  next  word, 
i.e.  two  bytes  later.  So  I  was  losing  data.  I  could  tell  I  was  having 
problems  because  I  was  comparing  it  to  a  known  file. 

24  bit  POSITION 
(bit  pos:  0..23) 


These  24  bits  hold  a  count  (in  16  bit  words)  of  user  clocks.  For 
example,  an  error  in  the  100th  serial  bit  of  a  serial  stream  would  have 
POSITION  count  of  6  indicating  an  error  presence  in  the  7th  16  bit 
word.  To  find  which  exact  bit  of  the  word  was  in  error,  one  would  look 
at  which  bit  in  the  XOR  word  was  a  logical  "1"  (mismatch).  In  this 
example,  it  would  be  the  4th  bit.  It  is  guaranteed  that  at  least  one 
event  will  occur  per  the  16,777,218  counts  of  the  24  bit  counter  to 
allow  a  software  implementation  of  a  larger  counter. 

( The  guaranteed  event  promised  to  detect  the  clock  rollover  is  that  the 
24  bit  counter  reads  2,  and  the  XOR  word  is  all  zeros.  So  all  I  did  to 
detect  clock  roll  over  was  to  see  if  one  the  current  clock  value  was 
less  than  the  previous.) 


bits  EVENT 
(bit  pos;  24. .31) 

Four  of  the  eight  bits  in  this  byte  are  modifiers  to  the  event  type. 

These  modifiers  are  orthogonal  in  that  one  event  could  have  several 
modifiers.  For  example,  an  event  could  have  errors  AND  be  a  used  as  a 
marker  AND  be  a  RESYNC  location.  The  other  four  bits  in  the  event  byte 
are  used  to  indicate  the  exact  bit  position  in  which  the  marker  event 
occurred.  THe  marker  position  bits  are  only  valid  when  the  Marker  event 
bit  is  set  to  1 .  The  modifiers  are  defined  as  follows: 


bO  (bit  pos:  24) 
bl  (bit  pos:  25) 
b2  (bit  pos:  26) 
b3  (bit  pos:  27) 
b4-b7(bit  pos:  28-31) 
word. 


M-Marker  (Rear  Panel) 

P-Pattem  Cycle  Marker  (Internal) 

R-Start  Resynchronization  Event 
S-Begin  Error  Squelch  Event 
Bit  position  of  the  Marker  within  this  event 


Marker:  This  is  an  internal  pole  once  every  second  (for  Edward's  Air 
Force  base).  We  didn't  worry  about  this  event. 

Pattern:  The  BitAlyzer  has  acquired  the  pattern,  and  is  now  synched  up. 
NOTE:  The  PATTERN  event  might  only  happen  once  in  a  file,  Usually  at 
the  beginning  of  a  file. 

Resync:  The  BitAlyzer25  has  reached  a  threshold  of  errors  where  it  no 
longer  is  able  to  hold  pattern,  and  now  is  in  the  process  of  trying  to 
acquire  pattern. 

Squelch:  To  many  errors  have  occurred  and  have  filled  the  BitAlyzer 
internal  buffer.  The  BitAlyzer  empties  the  buffer,  and  then  resumes 
detecting  and  writing  errors  to  the  file.  During  the  squelch,  errors 
are  NOT  recorded.  The  buffer  is  2  k  words,  and  is  controlled  by 
software. 


Okay,  now  here  comes  some  iffy  stuff.  When  I  asked  them  how  to  detect 
when  the  events  were  over,  they  told  be  the  following.  They  said  that 
for  the  Resync  or  Squelch  to  end,  the  ZERO  event  had  to  occur.  The  zero 
event  was  defined  to  be  that  the  XOR  word  was  all  zeros,  and  that  the 
EVENT  Byte  had  to  be  all  zeros.  I  used  this  in  my  code,  but  when  I 
would  compare  it  to  known  data  I  could  never  match  up  the  results.  So  I 
used  the  fact  that  they  said  that  no  errors  were  recorded  during  a 
Squelch  event,  and  some  trial  and  error  to  decide  for  myself  that  the 
way  to  detect  that  the  events  are  over  if  bits  24-27  are  all  zero,  i.e. 
the  marker  events  are  zero.  I  did  this  and  presto,  my  data  matched  the 
output  their  expensive  code  said  that  it  should.  I  still  have  a  small 
discrepancy  on  the  number  of  bits  lost  in  a  squelch  event,  but  the 
correct  number  of  squelch  events  is  calculated.  I'm  not  sure  if  the 
discrepancy  comes  from  the  fact  that  I  window  the  data  or  not.  It's  been 
6  months  since  I've  worked  on  the  problem,  and  my  memory  is  rusty  on  the 
details.  The  real  trick  to  figuring  out  the  event  structure  was  to 
compare  it  to  a  known  file. 

The  man  to  contact  at  Synthesis  Research  is  a  guy  called  Bert  Camer 
(650)  364-1853,  Bert_Camer@synthesysresearch.com.  Tell  him  that  you're 
working  with  Edward's  Air  Force  Base,  and  he'll  probably  be  helpful. 

*/ 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  "ba25read.h" 

#ifdefMATLAB 

#include  "D:\MATLABR1  l\extem\include\mex.h" 

#endif // ifdef  MATLAB 

/* 

*  Declare  local  routines. 

*/ 

void  ProcessErrorFile(  FILE  *fin,  FILE  *fdata,  FILE  *fout,  double*,  double*,  double*) ; 
/* 

How  the  ErrsPerByte[]  array  is  used  to  find  the  number  of  I's: 
ErrsPerByte[BinaryNumberIndex]=No.  of  I's  in  the  BinaryNumberIndex 
ErrsPerByte[5]— >ErrsPerByte[b0101]— >  two  I's  =  2; 

*/ 

unsigned  char  EiTsPerByte[]  =  { 

0,  1, 1,2,  1,2,  2,  3, 1,2,  2,  3, 2,  3,  3, 4, 

1.2,  2,  3,  2,  3,  3,  4,  2,  3,  3,  4,  3,  4,  4,  5, 

1.2. 2,  3,  2.  3.  3,  4, 2,  3,  3, 4,  3, 4, 4,5, 

2,  3,  3,  4,  3,  4,  4,  5,  3,  4,  4,  5,  4,  5,  5,  6, 

1.2,  2,  3,  2,  3,  3,  4,  2,  3,  3,  4,  3,  4,  4,5, 


2,  3, 3, 4,  3, 4, 4,  5,3,4, 4,  5, 4,  5,5,  6, 

2,  3,  3,  4,  3,  4, 4,  5,3,4,  4,  5, 4,  5,  5,  6, 

3,  4,  4,  5,  4,  5,  5,  6,  4,  5,  5,  6,  5,  6,  6,  7, 

1,2,  2,  3,  2,  3, 3,  4,  2,  3,  3, 4,  3, 4,  4,5, 

2,  3,  3,  4,  3, 4,  4,  5,3,4,  4,  5,  4,  5,  5,  6, 

2,  3,  3,  4,  3,4,  4,  5,3,4,  4,  5,  4,  5,  5,  6, 

3.4. 4.5.4,  5,5,6, 4,  5,  5,  6,  5,  6,  6,  7, 

2,  3,  3,  4,  3, 4,  4,  5,3,4,  4,  5, 4,  5,  5,  6, 

3. 4. 4,  5,  4,  5,  5,  6, 4,  5,  5,  6,  5,  6,  6,  7, 

3.4,  4,5,4,  5,5,6,  1,  5,  5,  6,  5,  6,  6,  7, 

4,  5,  5,  6,  5,  6,  6,  7,  5,  6,  6,  7,  6,  7,  7,  8}  ; 


/* 

Routine  Definitions 

*/ 

#ifndefMATLAB 
void  niain(  int  argc.  char  **  argv  ) 

{ 

FILE  *  errorFile  ; 

FILE  *  dataFile  ; 

FILE  *  outFile ; 

char  fileXfer[MAX_XFER_SIZE]; 

int  numBytes; 
int  numCycIes=0; 
int  i=0; 

int  actErrors=0; 
double  BER; 
double  ActualErrs; 
double  TotalBiis; 


/* 

Open  the  BA25  Error  File  and  process  the  events. 
*/ 

if(argv[l] ) 

{ 


errorFile  =  fopen(  argv[l],  "rb") ; 
dataFile  =  fopen(  argv[2],  "rb"); 
outFile  =  fopen(  argv[3],  "wb+"); 

while  ( !feof(dataFile) ) 

{ 

//  copy  the  data  file  to  the  ouput  file 

numBytes  =  fread(fileXfer,l,MAX_XFER_SIZE, dataFile); 

fwrite(fileXfer,  1  ,numBytes,outFile); 

numCycles++; 
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} 


rewind(dataFile); 

rewind(outFile); 

if  (  errorFile  &&  dataFile  &&  outFile ) 

ProcessErrorFile(  errorFile,  dataFile,  outFile,  &BER,  &ActualErrs, 

&TotalBits ); 

fclose(  errorFile ); 
fclose(  dataFile ); 
fclose(  outFile ); 

} 

} 

else 

printfC'Syntax:  procBERA  BERAfile  binarydatafile  outputfile\n\n"); 

} 

#endif 


*  Function  :  void  ProcessErrorFile(  FILE  *  fin,  FILE  *fdata,  FILE  *fout, 

double*  BER,  int* 


ActualErrs,  int*  ActualBits ) 


*  Purpose/Descr. 

*  Count  the  total  number  of  bits  processed  and  the  total  number  of 

*  errors  recorded.  Also  xor  the  errors  with  the  original  data  file 

*  to  produce  a  new  file  with  errors  in  the  positions  described  by  the 

*  error  file. 

*  Input: 

*  FILE  *  fin  -  pointer  to  the  error  event  file  to  be  processed 

*  FILE  *  fdata  --  pointer  to  the  original  error-free  data  file 

*  FILE  *  fout  -  pointer  to  a  replica  of  the  data  file.  This  file  will 

be  modified  with  the  bit  errors 

*  double  *  BER  -  double  pointer  to  variable  that  will  receive  the  actual  Bit  Error  Rate 

*  int  *  ActualErrs  -  int  pointer  to  the  variable  that  will  receive  the  actual  number  of  errors 

*  int  *  ActualBits  -  int  pointer  to  the  variable  that  will  receive  the  actual  number  of  bits 

processed 

*  Output:  void 

*  Return  Codes:  None 

He*************************************************************************/ 


void  ProcessErrorFile(  FILE  *fin,  FILE  *  fdata,  FILE  *fout, 
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double*  BER,  double*  ActualErrs,  double*  ActualBits ) 


unsigned  long 

LastPosition  =  Oxffffffff; 

unsigned  long 

PositionCounter; 

unsigned  short 

EventByte; 

unsigned  short 

Xor; 

unsigned  long 

Efi; 

Sample25 

ThisSample; 

double 

TotalBits=0; 

double 

TotalErrors=0; 

unsigned  long 

RollOvers=0; 

Word 

ThisWord; 

char 

fileCompl[MAX  XFER  SIZE]; 

char 

fileComp2[MAX_XFER_SIZE]; 

int 

numBytes; 

int 

numCycles=0; 

int 

i; 

rewind(fdata); 

rewind(fout); 

//  read  one  record  (16-bit  word)  from  the  error  file 
while(  fread(  &ThisSample,  sizeof(Sample25),  1,  fin ))  //  16-bit  words 

{ 

/* 

Extract  the  PositionCounter,  Event,  and  Xor  fields  from  the 
error  event. 

*/ 

PositionCounter  =  ThisSample.Clock; 

EventByte  =  HW25_EVENT_BYTE(  &ThisSample  ), 
Xor  =  ThisSample.Xor; 


/* 

Check  for  Rollovers  on  the  PostitionCounter 

The  guaranteed  event  promised  to  detect  the  clock  rollover  is  that  the 
24  bit  counter  reads  2,  and  the  XOR  word  is  all  zeros.  So  a  rollover 
occurred  whenever  the  position  counter  is  less  than  the  last  value  of 
the  position  counter. 

*/ 

if(  (PositionCounter  <  LastPosition)  &&  (LastPosition  !=  Oxffffffff) ) 
RollOvers-!“i-; 


/* 

Process  the  Errors 
*/ 

if  (Xor) 
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/* 

Compute  the  error  free  interval.  "Efi"  represents  the  error  free 
interval  as  the  number  of  error-free  words.  Words  are  16  bits  long. 

If  (PositionCounter  -  LastPosition)  equals  1 , 

then  there  are  no  error-free  words  in  between  the  previous  event 

and  the  current  word. 

*/ 

if(  LastPosition  ==  Oxffffffff ) 

//  if  LastPosition  is  Oxffffffff  then  this  is  the  first  error 
//  so  the  error-free  interval  is  equal  to  the  PositionCounter 
Efi  =  PositionCounter ; 

else  if(  /* (PositionCounter  >  LastPosition)  &&  */!RollOvers ) 

//  else  if  there  haven't  been  any  rollovers,  then  the  error-free 
//  interval  is  simply  the  difference  between  the  two  counters 
Efi  =  PositionCounter  -  LastPosition  -  IL; 

else 

//  otherwise  at  least  one  rollover  has  occured,  so  the  math  is 
//  a  little  different 

Efi  =  RollOvers*ROLL_PERIOD_25  -  (LastPosition  - 

PositionCounter)  -  IL; 


//  only  keep  track  of  how  many  rollovers  occur  between  error  events 
//  so  reset  this  counter  every  time  an  error  is  detected. 

RollOvers  =  0; 


if  (!feof(fdata)) 

{ 

//  Seek  to  an  of  offset=Efi  from  the  current  position  in  the 
//  data  file.  That  is  where  the  next  error(s)  occurs, 
if  (fseek(fdata,Efi*sizeof(Word),SEEK_CUR)  “0) 

{ 

ftell(fdata); 

//  Should  read  2  bytes 

numBytes  =  fread(&ThisWord,l,sizeof(Word),fdata); 

if  (numBytes) 

{ 

//  Flip  the  errored  bits 
ThisWord.byteO  ((Xor  »  8)  &  Oxff); 
ThisWord.bytel  '^=  (Xor  &  Oxff); 

//  Go  to  the  same  position  in  the  ouput  file  (received 

data) 


fseek(fout,Efi*sizeof(Word),SEEK_CUR); 
//  Write  in  the  data  with  errors 
fwrite(&This  Word,  1  ,numBytes,fout); 
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error  count 


//  Add  the  number  of  errors  in  this  word  to  total 


if(numBytes>=  1 ) 

TotalErrors  +=  EiTsPerByte[  (Xor»8)  & 

Oxff]; 

if(numBytes>=2) 

TotalErrors  +=  ErrsPerByte[  (Xor&Oxff)  ]; 

/* 

Add  the  number  of  Error-Free  Words  (times  16 

bits/word) 

plus  the  number  of  bits  in  the  XOR  field  (16  bits) 
to  the  total  number  of  bits  processed. 

*/ 

TotalBits  +=  (Efi  *  16); 

TotalBits  +=  8*numBytes; 

} 

} 

else 

//  else  fseek  couldn't  locate  this  offset  (EOF) 
break; 

} 

else 

//  if  data  file  is  shorter  than  the  error  file,  then  ignore  the  rest 
//  of  the  error  file, 
break; 

} 

LastPosition  =  PositionCounter; 

} 

#ifdef  DEBUG 
/* 

Compare  the  output  file  to  input  file  to  determine  the  actual  number  of  errors 
*/ 

rewind(fout); 

rewind(fdata); 

*ActualErrs=0; 

*ActualBits=0; 

while(!feof(fdata)  &&  !feof(fout)) 

{ 

//read  in  a  big  chunk  of  the  file  each  time 
ftell(fout); 

numBytes  =  fread(fileCompl,l,MAX_XFER_SIZE,fout); 
fread(fileComp2, 1  ,MAX_XFER_SIZE,fdata); 
for(i=0;  i<numBytes;  i++) 

{ //  compare  input  and  output  data  byte-by-byte 

*ActualErrs  +=  ErrsPerByte[(fileCompl[i]  ^  fileComp2[i])&0xff]; 
*ActualBits  +=  8; 


129 


} 


} 

#else 

*ActualErrs  =  TotalErrors; 
fseek(fout,OL,SEEK_END); 

*ActualBits  =  ftell(fout)*8; 

#endif 

*BER  =  (*ActualErrs)/(*ActualBits); 

/*  printf(  "Total  Pits  Processed  =  %.01f\n",(float)(*ActualBits)); 
#ifdef  DEBUG 

printf(  "Total  Errors  Counted  =  %.01f\n",TotalErrors); 

#endif 

printf(  "Total  Actual  Errors  =  %.01f\n",(float)(*ActualErrs)); 
printf(  "Bit  Error  Rate  =  %6.4e\n",*BER); 

*/ 

} 


#ifdefMATLAB 

//  matlab  gateway  function  makes  it  possible  to  call  this  function  from  matlab. 

//  to  use  this  function  in  matlab,  compile  it  on  the  matlab  command  line  as  follows: 
//  1 .  change  the  working  directory  path  to  D:\stats\simulink  simulation\ 

//  »  cd  'D:\stats\simulink  simulation\' 

//  2.  compile  this  program  as  a  mex  file: 

//  »  mex  -DMATLAB  -output  procBERA  procBERA.c 

//  3.  use  this  function  in  a  MATLAB  script  as  a  normal  function  call 

//  »procBERA('d:\stats\flight2bera\F02f_02_M12_O28.bin', 

'd:\stats\SDS00 1  c.dat',  'd:\stats\rcvddata.bin',BER,actualErrs,actualBits) 

void  mexFunction(int  nlhs,  mxArray  *plhs[],int  nrhs,  const  mxArray  *prhs[]) 

{ 

int  buflen; 
char*  errfile; 
char*  datafile; 
char*  outfile; 

FILE*  errorFile; 

FILE*  dataFile; 

FILE*  outFile; 

char  fileXfer[MAX_XFER_SIZE]; 

int  numBytes; 

int  numCycles=0; 

double*  BER; 

double*  ActualErrs; 

double*  TotalBits; 

/*  Check  for  proper  number  of  arguments.  */ 
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if(nrhs!=3) 

{ 

mexErrMsgTxt("Three  inputs  required."); 

} 

else  if(nlhs>3)  { 

mexErrMsgTxt("Too  many  output  arguments"); 

} 

/*  The  inputs  must  be  strings  */ 

if(  !mxIsChar(prhs[0])  1|  !mxIsChar(fvlis[l])  ||  !mxIsChar(prhs[2])) 

{ 

mexErrMsgTxt("Inputs  must  be  strings."); 

} 

/*  Get  the  length  of  the  input  string.  */ 

buflen  =  (mxGetM(prhsp])  *  mxGetN(prhs[0]))  +  1; 

/*  Allocate  memory  for  input  and  output  strings.  */ 
errfile=mxCalloc(buflen,  sizeof(char)); 

/*  Assign  pointers  to  each  input  and  output.  */ 
mxGetString(prhs[0],errfile, buflen); 

/*  Get  the  length  of  the  input  string.  */ 

buflen  =  (mxGetM(prhsn])  *  mxGetN(prhs[l]))  +  1; 

/*  Allocate  memory  for  input  and  output  strings.  */ 
datafile=mxCalloc(buflen,  sizeof(char)); 

/*  Assign  pointers  to  each  input  and  output.  */ 
mxGetString(prhs[  1  ],datafile,buflen); 

/*  Get  the  length  of  the  input  string.  */ 

buflen  =  (mxGetM(prhsp])  *  mxGetN(prhs[2]))  +1; 

/*  Allocate  memory  for  input  and  output  strings.  */ 
outfile=mxCalloc(buflen,  sizeof(chru)); 

/*  Assign  pointers  to  each  input  and  output.  */ 
mxGetString(prhs[2],outfile, buflen); 

/*  Create  matrix  for  the  return  argument.  */ 
plhs[0]  =  mxCreateDoubleMatrix(l,l,  mxREAL); 
plhs[l]  =  mxCreateDoubleMatrix(l,l,  mxREAL); 
plhs[2]  =  mxCreateDoubleMatrix(l,l,  mxREAL); 

BER  =  mxGetPr(plhs[0]); 

ActualErrs  =  mxGetPr(plhs[l]); 

TotalBits  =  mxGetPr(plhs[2]); 

errorFile  =  fopen(  errfile,  "rb") ; 
dataFile  =  fopen(  datafile,  "rb"); 
outFile  =  fopen(  outfile,  "wb+"); 


//  Copy  the  data  file  to  the  output  file 


while  ( !feof(dataFile) ) 

{ 

numBytes  =  fread(fileXfer,l,MAX_XFER_SIZE,dataFile); 

fwrite(fileXfer,  1 , numBytes, outFile); 

numCycles++; 

} 

rewind(dataFile); 

rewind(outFile); 

if  (  errorFile  &&  dataFile  &&  outFile ) 

ProcessEiTorFile(  errorFile,  dataFile,  outFile,  BER,  ActualErrs,  TotalBits  ); 
fclose(  errorFile ); 
fclose(  dataFile ); 
fclose(  outFile ); 

} 

else 

printf(" Syntax:  procBERA\(BERAfile,  binarydatafile,  outputfile\)\n\n"); 


#endif 


function  numErrs=errsPerByte (byte) 
if  (nargin  •'-=  1 ) 

error (' ErrsPerByte  musr  have  one  argument:  ErrsPerEyte ( Byte ) \n ' ) 
else 

errors^ [ .  .  . 

0,  1/  If  2,  If  2f  2f  3f  If  2,  2,  3f  2f  3f  3,  4f... 

1,  2f  2f  3,  2f  3,  3f  4f  2,  3,  3,  4,  3,  4f  4f  5f  .  .  . 

1,  2,  2f  3,  2f  3,  3,  4,  2,  3f  3f  4,  3,  4,  4,  5,... 

2,  3,  3f  4f  3,  4,  4f  5,  3,  4,  4f  5,  4,  5f  5,  6f... 

1,  2f  2f  3f  2,  3f  3f  4,  2,  3f  3,  4f  3,  4,  4,  5,... 

2f  3,  3f  4,  3f  4,  4,  5,  3,  4f  4,  5f  4,  5,  5,  6,... 

2,  3f  3f  4f  3,  4f  4f  5,  3f  4,  4,  5,  4,  5f  5,  6f  .  .  . 

3,  4,  4f  5f  4,  5f  5f  6,  4,  5,  5,  6,  5,  6,  6,  7,... 

1,  2f  2f  3,  2f  3,  3,  4,  2f  3,  3f  4,  3,  4,  4,  5,... 

2f  3f  3,  4f  3f  4f  4,  5f  3,  4f  A ,  5/  4,  5/  5f  6,... 

2f  3,  3,  4,  3,  4f  4f  5,  3f  4,  4f  5f  4,  5,  5,  6,... 

3f  4,  4f  5,  4f  5f  5,  6,  4,  5,  5,  6,  5,  6,  6,  7,... 

2f  3,  3,  4,  3f  4,  4f  5,  3,  4,  4f  5f  4,  5f  5,  6,... 

3f  4,  4f  5f  4f  5f  5f  6,  4f  5,  5,  6,  5,  6,  6,  7,... 

3,  4f  4f  5f  4f  5,  5f  6,  4,  5,  5,  6,  5,  6,  6,  7,... 

4,  5,  5,  6,  5,  6,  6,  7,  5,  S,  6,  7,  6,  7,  7,  8]; 


numErrs=errors (byte  +  1)  ; 

end 

function  [BERf  numErrs , numBits ] =checkErrorRate (infilOf  outfiie) 
?;checl:ErrorRare 

if (nargin  2 ) 

error ( ' checkError Rate  needs  2  arguments :  checkErrorRate (inril 
curfile) \n' ) ; 
else 

fin  =  f open ( inf ile, ' rb ’) ; 
fout  =  fopen (outfiie, ’ rb ') ; 
numErrs^O ; 

if({fin-=0)  &  (fout~=0)) 

while  (feof(fin)  ==  0) 

[input, numBytes] =fread ( fin, hex2dec ( ' iuOOUUO ’ )  ,  ' uchar ' ) ; 
output==fread(fout,  numBytes,  ’uchar'  )  ; 

numErrs  =  numErrs  + 

ErrsPerByte (bitand (bitxor (input { 1) , output { 1) ) , hex2dec ( ’ f f ' ) ) ) ; 

1  end 

end 

end 

end 
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Data  Statistics  Appendix 


The  following  pages  contain  the  statistical  analysis  for  the  probabilities  and  entropies  for 
each  of  the  data  sets. 
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Symbol  number  H=2.64 
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Symbol  number  H=2.8 
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Symbol  number  H=2.52 
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Symbol  number  H=  7,05 
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Symbol  number  H=  5.94 
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Symbol  number  H=  6.56 
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Symbol  number  H=  5.94 


Frequency  information  for  2 
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Symbol  number  H=  5.89 
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Symbol  number  H=  5.77 
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Symbol  number  H=  4.62 
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Symbol  number  H=  4.87 
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Symbol  number  H=3.71 
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Symbor  number  H=  3.72 


